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SECTION  I 


INTRODUCTION  AND  SUMMARY 


The  objective  of  this  effort  was  to  document  the  performance  capabilities 
of  an  Infrared  Variable  Acuity  Remote  Viewing  System  (VARVS)  and  to  optimize 
its  parameters.  In  addition,  the  feasibility  of  adding  a common  aperture 
narrow  f ield-of-view,  high  resolution  optical  system  to  the  wide  angle  VARVS 
was  to  be  determined. 

This  turned  out  to  be  a far  more  ambitious  program  than  we  originally 
anticipated.  An  entirely  new  observer  dynamic  model  had  to  be  developed 
because  no  existing  model  could  properly  handle  the  large  f ield-of-view  of 
the  VARVS  system.  The  next  difficulty  arose  when  we  decided  that  the  only 
meaningful  way  to  present  VARVS  capabilities  was  on  a comparative  basis  with 
conventional  linear  sensors.  This  required  that  we  develop  a model  for  the 
linear  sensor  along  lines  identical  to  the  VARVS,  so  common  inputs/outputs, 
etc.  would  exist.  This  in  turn  required  that  we  validate  both  models 
experimentally.  When  we  considered  that  a minimum  of  three  f ields-of-view 
are  required  to  validate  the  linear  sensor  model  and  the  one  f ield-of-view 
of  the  VARVS,  the  experimental  effort  expanded  by  a factor  of  four.  Fortu- 
nately much  of  the  extra  experimental  work  was  able  to  make  use  of  test 
facilities  developed  for  validation  of  electro-optic  sensor  models  that  were 
part  of  an  on-going  Independent  Research  and  Development  (IRAD)  effort.  As 
a result,  the  contract  objectives  were  met  within  available  funding  levels, 
but  a slip  in  schedule  necessarily  occurred. 

The  results  of  this  study  were  somewhat  negative.  The  Infrared  VARVS 
target  detection  performance  was  found  to  be  at  best  only  80%  as  effective 
as  the  conventional  linear  sensor.  The  limiting  factor  is  the  VARVS  slew 
time.  It  becomes  dominant  because  its  optical  axis  must  be  slewed  for  each 
cue  before  the  acuity  is  sufficient  for  target  detection.  Conversely,  the 
linear  sensor's  uniform  resolution  supports  search  and  detection  throughout 
its  f ield-of-view.  The  major  portion  of  visual  search  can  then  be  accomplished 
by  visual  glimpses  with  sensor  slewing  only  occassionally  dictated  by  lateral 
coverage  requirements.  VARVS  performance  can  be  improved  either  by  reducing 
its  slew  time  and  including  control  by  eye  position  or  by  increasing  its 
region  of  high  acuity.  The  latter  approach  results  in  continual  improvement 
until  the  entire  field  is  uniform  where  it  then  becomes  a linear  sensor  with 
about  20°  f ield-of-view. 


An  obvious  question  at  this  point  is,  "Why  use  the  VARVS  when  a linear 
sensor  works  better?".  For  target  acquisition  this  is  clearly  the  case.  It 
was  shown  both  experimentally  and  theoretically  in  this  effort.  Why  then 
is  VARVS  display  so  appealing  when  viewing  real  world  scenes?  This  is 
especially  true  when  video  tapes  made  in  flight  are  viewed  but  also  exists  to 
a lesser  degree  when  viewing  from  a laboratory  rooftop.  The  answers  probably 
lie  in  the  cognitive  factors  involved  in  the  visual  process  that  are  not 
defined  or  well  understood.  These  include  peripheral  cueing  due  to  motion  or 
glints.  Others  may  be  the  instantaneous  recognition  of  relative  vehicle, 
sensor,  real  world  orientation  that  is  possible  with  a wide  field  display. 

This  may  be  especially  valuable  when  the  vehicle  is  maneuvering  and/or  flight 
control  and  navigation  must  be  performed  simultaneously  with  target  acquisition. 
We  believe  therefore,  that  neither  the  theory  or  experiment  show  the  true  value 
of  the  VARVS.  It  is  becoming  increasingly  clear  that  only  a flight  evaluation 
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One  significant  technical  output  of  this  contract  is  the  new  observer 
dynamic  model.  This  model  synthesizes  human  vision  much  more  accurately  than 
existing  models.  The  primary  difference  is  shown  in  Figure  1.  This  new  model 
when  applied  to  conventional  sensors  can  be  used  to  properly  evaluate  advanced 
wide  angle  display  concepts  whereas  the  existing  models  are  only  verified 
for  10°  to  20°  display  angles. 

The  feasibility  of  a dual  field  VAE.VS  could  not  be  positively  established 
during  this  study.  However,  results  were  encouraging  enough  that  we  Intend 
to  continue  this  effort  independently. 
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SECTION  II 
APPROACH 


The  study  was  divided  into  two  phases.  The  primary  objective  of  Phase 
I was  to  develop  an  experimentally  validated  analytic  model  of  the  VARVS 
concept  and  to  use  this  model  to  predict  target  acquisition  performance  of 
an  infrared  VARVS.  The  secondary  objective  was  to  determine  the  amount  of 
variation  that  could  be  made  in  the  VARVS  design  parameters  that  would 
yield  improved  performance.  In  Phase  II,  the  feasibility  of  the  improved 
VARVS  was  to  be  established.  As  the  study  progressed,  it  became  clear  that 
the  VARVS  model  development,  its  experimental  validation,  and  its  application 
to  the  IR  VARVS  would  occupy  most  of  the  effort.  This  occurred  because 
all  existing  sensor  dynamic  models  were  not  useable  with  the  VARVS  because 
they  are  developed  around  linear  sensor  and  display  geometries.  Furthermore, 
to  avoid  the  pitfalls  of  predicting  absolute  performance  of  the  VARVS,  we 
elected  to  show  performance  relative  to  a conventional  linear  sensor.  This 
required  a parallel  development  of  a linear  sensor  model  based  on  the  same 
philosophy  as  the  VARVS  model  so  that  it  would  employ  the  same  inputs  and 
outputs.  This  linear  model  has  been  developed  and  experimentally  validated 
under  a related  IRAD  effort,  see  Reference  1,  and  could  be  used  for  this 
purpose. 

The  actual  approach  used  for  Phase  I is  shown  in  Figure  2.  As  shown  on 
the  figure,  the  study  began  with  an  effort  to  improve  the  VARVS  static  model, 
using  as  a base,  the  experience  gained  during  the  development  of  the  Office 
of  Naval  Research  (ONR)  brassboard  system.  More  specifically,  the  acuity 
function  and  optical  quality  inputs  to  the  VARVS  model  were  updated  to 
reflect  actual  performance  of  the  ONR  brassboard  equipment  that  was  to  be 
used  for  the  validation  experiment.  The  major  effort  was  directed  to  the 
inclusion  of  the  human  observer  performance  to  create  the  dynamic  model  for 
the  VARVS.  This  effort  is  described  in  detail  in  Section  III.  In  parallel 
with  this  task  but  funded  by  a related  IRAD,  was  the  development  of  a con- 
ventional sensor  model  along  identical  lines.  The  design  of  the  experimental 
equipment  was  also  accomplished  at  this  point  in  the  study  which  consisted 
of  a treadmill  and  drive  system  to  simulate  relative  forward  vehicle  motion. 
Details  of  this  development  are  also  discussed  in  Section  III.  The  next 
step,  the  analytic  model  validation,  was  accomplished  by  inputting  a set  of 
flight  conditions  into  the  model  and  comparing  the  predicted  performance 
with  actual  results  from  the  experiment  for  the  same  simulated  flight 
conditions.  Correlation  of  theory  and  experiment  led  to  several  iterations 
of  the  theoretical  model  before  agreement  was  achieved. 

For  the  parametric  analyses  discussed  in  Section  IV,  both  conventional 
and  VARVS  models  were  shifted  to  the  8-14y  infrared  spectrum  so  that  perform- 
ance on  an  Infrared  VARVS  (IR  VARVS)  could  be  compared  directly  to  a FLIR 
with  conventional  optics.  Both  IR  VARVS  and  FLIR  performance  were  computed 
as  functions  of  the  mission  parameters  of  altitude,  velocity,  minimum 
acquisition  range  (for  weapon  delivery  considerations),  and  lateral  coverage. 
In  addition,  the  effect  of  the  IR  VARVS  slew  time  and  FLIR  f ield-of-view 
performance  were  determined.  Slew  time  was  included  because  it  has  the 
greatest  effect  of  any  parameter  on  the  VARVS  performance.  A separate  off- 
shoot of  this  parametric  analysis  was  a study  of  the  internal  structure  of 
the  VARVS  model  to  see  if  an  alternate  distortion  function  would  yield  better 
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Figure  2 Phase  I approach 


performance.  This  study  Is  discussed  In  Section  V while  details  of  the  IR 
VARVS/FLIR  parametric  study  are  presented  In  Section  IV. 

Rationale'  used  for  performance  comparisons  of  the  IR  VARVS  to  FLIR  Is 
shown  in  Figure  3.  First,  the  probability  of  target  detection  was  computed 
and  plotted  as  a function  of  offset  for  all  parameter  combinations.  Next, 
the  lateral  coverage  where  the  probability  of  detection  (Pd)  la  90X  was 
read  from  each  curve.  From  this  point,  the  approach  depends  upon  the  type 
of  mission.  For  a reconnaissance  mission,  the  90Z  lateral  coverage  was 
plotted  as  a function  of  depression  angle  for  each  parameter  set.  For 
weapon  delivery,  the  90Z  P<j  lateral  coverage  was  plotted  against  the 
minimum  range  boundary  of  the  downrange  search  area.  For  reconnaissance, 
the  maximum  lateral  coverage  was  selected  from  each  curve.  A ratio  was 
then  computed  of  IR  VARVS  and  10*,  20*,  and  40*  fleld-of-vlew  FLIR  lateral 
coverage  for  Identical  parameter  sets.  These  data  were  plotted  against  the 
remaining  parameters  which  are  VARVS  slew  time,  clutter  factor,  altitude, 
and  velocity.  This  technique  of  presentation  was  selected  to  show  Instantly 
the  relative  performance  of  sensor  types  and  sensitivities  to  design  and 
mission  parameters.  Similar  curves  were  plotted  for  weapons  delivery 
except  In  this  case  the  product  of  lateral  coverage  at  maximum  acquisition 
range  times  maximum  acquisition  range  was  plotted  as  the  dependent  variable. 
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Originally,  we  Intended  to  plot  the  ratio  of  IR  VARVS  to  FLIR  maximum 
acquisition  ranges  for  a specific  lateral  coverage  as  a figure  of  merit. 

This  proved  to  be  uninformative  because  the  ratios  were  simply  those  Imposed 
by  the  resolution  limits  of  the  sensors  until  the  lateral  coverage  caused 
a dynamic  search  limitation  which  then  caused  the  ratio  to  drop  abruptly 
to  zero.  For  these  reasons,  the  range  lateral  coverage  product  was 
selected  as  an  independent  variable  for  weapohs  delivery  plots. 

In  Phase  II  the  feasibility  of  an  improved  IR  VARVS  is  considered  by 
showing  how  an  on-axis  target  recognition  capability  can  be  added  to  the 
VARVS.  Details  are  presented  In  Section  VI. 
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SECTION  III 


ANALYTIC  MODEL  DEVELOPMENT  AND  VALIDATION 


In  this  section  the  development  of  the  analytic  models  Is  described 
and  their  experimental  validation  la  presented.  Since  validation  will  be 
accomplished  with  the  ONR  brassboard  equipment,  the  model  Is  first  developed 
for  visual  spectrum  operation.  Its  conversion  to  the  Infrared  la  described 
In  Section  IV. 

1.  STATIC  MODEL  DEVELOPMENT  «' 

A static  model  for  the  VARVS  has  been  developed  in  previous  efforts 
for  the  infrared  spectrum.  Reference  (2)  and  (3).  This  model  was  developed 
using  the  rationale"  developed  by  the  Army  Night  Vision  Laboratory  (NVL) 
which  is  undoubtedly  the  best  approach  currently  available.  The  adaptation 
of  the  infrared  NVL  model  (Reference  (4))  to  the  visual  spectrum  VARVS  is 
described  in  Appendix  A.  A simplified  flow  diagram  of  the  visual  VARVS 
model  Is  described  In  Figure  4.  The  model  uses  a slgnal-to-nolse  approach 
and  is  distinguished  by  the  use  of  a matched  filter  to  simulate  the  observer's 
pattern  recognition  capability  which  results  in  an  improved  threshold  of 
detection. 

The  performance  predictions  of  the  static  model  do  not  take  into  account 
the  time  constraints  under  which  the  observer  must  function.  Therefore, 
the  following  dynamic  model  was  developed  to  include  these  effects. 

2.  DYNAMIC  MODEL  DEVELOPMENT 

By  far  the  most  significant  analytical  effort  was  the  development  of 
the  VARVS  dynamic  model  and  its  companion  linear  sensor  model.  Development 
of  new  models  was  required  because  existing  models  do  not  include  the 
effects  of  peripheral  vision  which  is  obviously  very  important  when  modeling 
wide  field-of-view  displays  of  the  VARVS  type.  Conventional  linear  models 
then  become  unacceptable  for  the  linear  sensors  in  this  study  because  they 
require  empirical  inputs  such  as  clutter  factor  and  which  are  not  needed  for 
the  VARVS  model.  Fortunately  the  linear  sensor  model  was  being  developed 
under  IRAD  funding.  That  effort  is  described  in  Appendix  B.  The  VARVS 
dynamic  model  is  described  below. 

The  VARVS  system  provides  the  user  with  an  extremely  wide  field-of-view. 
The  resolution  across  this  field-of-view  is  primarily  a function  of  the 
foveal  lens  characteristic  distortlQn.  It  is  up  to  the  user  to  decide  how 
to  optimally  use  this  display  for  target  acquisition.  The  acquisition  process 
using  the  foveal  lens  is  a very  adaptive  process.  The  user,  as  will  be 
investigated  during  the  course  of  this  section,  must  decide  which  portion 
of  the  footprint  will  be  used  for  search,  and  how  he  will  spend  the  time  he 
has  available. 

The  following  is  a description  of  the  VARVS  target  acquisition  model. 

The  first  part  is  concerned  with  a general  description  of  the  VARVS  ground 
coverage  and  how  its  characteristics  affect  the  way  it  is  used.  The 
second  part  describes  the  exponential  equations  used  to  define  the  condlt- 
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Figure  4 Static  performance  model 


lonal  probabilities.  In  the  third  part  with  the  help  of  a flow  diagram, 
a general  description  of  the  total  process  is  given.  The  model  described 
below  is  the  result  of  numerous  iterations,  many  of  them  inspired  by  lack 
of  agreement  with  experimental  data. 

a.  General  Description 

The  extent  of  the  VARVS  ground  coverage  and  the  quality  of  the  coverage 
is  defined  in  terms  of  a sensor  footprint.  This  footprint  is  calculated  by 
partitioning  the  ground  plane  into  uniform  patches  (AA) . At  each  AA. 
sensor  performance  is  calculated.  The  performance  is  expressed  in  terms 
of  spatial  resolution  at  the  focal  plane  of  the  foveal  lens.  This  resolution 
is  calculated  as  a function  of  such  parameters  as  altitude,  velocity,  angle 
from  optical  axis,  atmospheric  conditions  and  internal  sensor/display  char- 
acteristics. These  calculations  are  accomplished  using  an  adaptation  of  the 
NVL  model  as  described  in  Reference  (4).  Along  with  the  resolution  calcula- 
tions made  at  each  AA,  the  size  of  the  cue  or  target  as  it  appears  on  the 
non-linear  lens  focal  plane  is  calculated.  The  product  of  these  two  terms, 
resolution  in  line  pairs/mm  and  cue  size  in  mm,  defines  the  number  of  line 
pairs  across  the  target.  Line  pairs  across  the  target  will  be  used  in 
conjunction  with  exponential  equations  similar  to  Johnson's  Detection  Criteria, 
Reference  (5)  to  predict  the  probability  of  detection  on  the  condition  that 
the  observer  is  looking  at  the  target. 

A typical  VARVS  sensor  footprint  defining  the  probability  of  detecting 
a 60  ft.  cue  is  shown  in  Figure  5.  As  is  expected  from  the  hature  of  the 
lens  acuity  function,  the  highest  probabilities  of  detection  are  present  on 
or  near  the  optical  axis,  and  the  quality  falls  off  rapidly  with  increasing 
angle  from  the  lens  optical  axis.  A somewhat  slmplier  description  of  the 
footprint  can  be  made  by  enclosing  regions  of  equal  probability  as  shown  in 
Figure  6 where  the  footprint  is  cut  up  into  a set  of  contours  of  equal 
probability. 

When  the  VARVS  is  used  for  target  acquisition,  the  operator  must  pick 
which  portion  of  the  sensor  footprint  he  will  use  for  the  search.  There  are 
two  main  forces  which  dictate  this  choice.  As  the  portion  of  the  footprint 
used  for  target  acquisition  increases,  it  becomes  easier  to  cover  the  total 
search  area.  At  the  same  time,  the  probability  of  detecting  the  target 
determined  by  imaging  over  the  portion  of  the  footprint  used,  decreases. 

The  target  is  less  likely  to  be  located  in  the  high  resolution  portion  of 
the  footprint  which  accounts  for  the  lower  average  probability. 

The  time  available  (t  ) to  find  a target  is  related  to  the  time  a 
target  stays  within  the  probability  contour  which  is  the  distance  (X  ) along 
the  flight  direction  of  the  chosen  search  footprint  contour  and  divined  by 
the  velocity  (Vx)  of  the  aircraft,  viz: 

X 

t - -E  (1) 
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Figure  6 Typical  probability  footprint 
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Detection  Contours 

Flight  Direction  ow-n»-t 

Figure  6 Probability  contours  of  VAR  VS  footprint 


In  order  to  conduct  an  effective  search  the  operator  must  cover  as 
much  of  the  search  area  as  possible.  This  may  be  accomplished  by  slewing 
the  sensor  from  side  to  side  in  a systematic  fashion,  taking  care  not  to 
spend  too  much  or  too  little  time  at  each  slew  position.  Both  slewing 
the  sensor  and  searching  the  display  take  time.  Specifically,  let  ts^  be 
defined  as  the  average  time  it  takes  to  slew  the  sensor  and  let  it  settle  to 
its  steady  state,  and  tg  be  the  time  related  to  a glimpse  at  the  display. 

The  time  per  glimpse  has  been  measured  to  be  about  0.25  to  0.3  seconds. 

The  operator  can  optimize  his  search  process  by  properly  adjusting 
the  allocation  of  time  to  sensor  slewing  and  the  time  spent  searching  the 
display  between  slews.  The  total  time  (t')  required  for  search  is 
calculated  in  Equation  (2). 


t' 

a 


«n 


V + t8l) 


Si 


(2) 


where  n defines  the  number  of  glimpses  per  slew  and  n ^ is  the  number  of 
slews  tftat  will  be  made  during  the  time  the  sensor  platform  travels  a 
distance  Xp.  Under  the  assumptions  of  a nominally  constant  depression 
angle 
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Thus: 
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The  upper  bound  on  t'  further  defines  the  restrictions  of  choices  of 

n and  n , chosen  for  optimal  results, 
g si 


b.  Probability  Calculations 


The  total  target  acquisition  probability  (P^,)  can  be  modeled  as  the 
product  of  four  probabilities,  viz: 


(5) 


where  P is  the  probability  of  detecting  the  true  cue  given  the  cue  is 
within  the  FOV  of  the  sensor,  Pd  is  the  probability  that  the  observer  will 
detect  the  target  if  he  is  looking  at  the  cue,  and  Pr  is  the  probability 
that  the  observer  will  recognize  the  target  if  he  is  looking  at  the  target. 
Finally,  P8  is  the  probability  that  the  true  cue  falls  within  the  chosen 
probability  contour  during  the  search  process.  The  probabilities  Pc,  Pd 
and  Pr  are  functions  of  line  pairs  across  the  target  or  cue. 


Now  we  will  define  the  probability  of  detecting  the  cue,  given  the 
observer  is  looking  at  the  cue  (P').  The  prime  indicates  that  this 
probability  does  not  take  into  account  multiple  glimpse  and  multiple  cue 
possibilities.  This  conditional  probability  of  cue  detection  (P')  is 
modeled  using  the  form, 
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(6) 


P'  ■ 1 - exp  [-(a  Nc  + b) 


2] 


where  Nc  Is  the  number  of  resolution  elements  across  the  cue  averaged  over 
the  portion  of  the  footprint  used.  Nc  has  units  of  line  pairs.  One 
resolution  element  across  the  cue  is  generally  accepted  as  the  threshold 
for  902  probability  of  target  detection.  It  has  also  been  found  that  the 
probability  of  detection  goes  to  zero  somewhere  In  the  range  of  1/7  to  1/8 
line  pairs  across  the  cue.  The  constants  a and  b are  evaluated  such  that 


?'  - 0.9  for  N - 1.0 
c c 

P'  - 0.0  for  N - 0.13 
c c 


(7) 


Using  Equation  (7)  to  define  a and  b In  Equation  (6)  results  In 

P'  - 1.  - exp  [-(1.75  Nc  - . 23) 2]  (8) 

Now  the  conditional  probability  of  target  detection  (P<j)  will  be 
defined.  It  Is  similar  to  P£  except  that  Instead  of  the  number  of 
resolution  elements  across  the  cue,  the  number  of  resolution  of  elements 
across  the  target  (N()  Is  used.  Thus, 

Pd  - 1 - exp  (-(1.75  Nt  - .23)2)  (9) 

The  third  probability  of  the  target  recognition  from  Bailey  (6)  Is: 


1 - exp  (-(- 


- D2) 


(10) 


A P -0.90  corresponds  to  five  line  pairs  across  the  target. 


If  a large  portion  of  the  display  Is  used  for  search,  the  observer  would 
probably  have  very  little  difficulty  in  covering  the  total  search  arec  but 
the  average  value  of  Nc  would  be  low.  At  the  other  extreme,  if  only  the 
area  of  the  display  near  the  optical  axis  Is  used  for  search,  the  observer 
would  be  hard  pressed  to  adequately  cover  the  total  search  area  but  the 
average  value  of  Nc  calculated  would  be  larger  than  In  the  first  case. 

The  effects  of  multiple  glimpses  and  multiple  targets  was  studied  in 
detail  before  deciding  on  the  type  of  model  to  be  used.  The  performance  of 
the  observer  can  best  be  analyzed  by  answering  the  question,  "What  Is  the 
probability  P(K)  of  finding  exactly  K targets  In  N glimpses?".  This 
probability  is  usually  associated  with  Binomial  Distribution  and  Is  defined 
by: 


,N.  K „ vN-K 
P(K)  - (R)  p (1  - p)  , 
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In  the  multiple  cue  acquisition  case,  the  probability  of  finding  a 
cue  is  a function  of  the  number  of  unfound  cues  in  the  display;  therefore, 
as  the  observer  detects  cues,  the  probability  of  finding  another  cue 
changes.  The  equation  which  relates  the  probability  of  finding  a cue  per 
glimpse  for  one  cue  contained  within  the  display  (P£)  to  the  probability 
(P(NT))  of  finding  a cue  where  multiple  independent  (NT)  cues  exist  is: 


P(NT)  - 1 - (1  - P')NT  (12) 

As  the  observer  detects  and  responds  to  possible  cues,  P(NT)  decreases 
because  the  number  of  unfound  cues  is  decreasing.  Due  to  this  non-stationary 
characteristic  of  P(NT),  the  expected  number  of  cues  found  during  K glimpses 
E(NT)  cannot  be  calculated  using  compact  methods  such  as  Bernoulli  trials 
(a  form  of  the  Binomial  expansion) . Therefore  in  the  next  paragraphs  the 
probability  density  functions  to  cover  this  case  will  be  developed. 


The  probability  density  function  must  define  the  probability  of 
finding  K cues  in  nn  glimpses  for  constant  n^  and  for  all  K such  that  K is 
an  Integer  where 


g 


g 


0 < K < n 
- - g 

Then  one  must  determine  all  possible  ways  to  detect  K cues  in  ng  glimpses 
(all  permutations),  find  the  probability  associated  with  each  permutation, 
and  sum  these  probabilities.  This  process  is  simplified  by  using  a digital 
computer  and  its  memory  capabilities. 

Consider  the  state  diagram  in  Figure  7.  The  circles  represents  specific 
states  as  defined  by  the  two  axes,  number  of  glimpses  made  (trials)  and  cues 
detected  (successful  results).  Each  state  has  associated  with  it  the 
probability  of  being  in  that  state.  The  probability  of  traveling  from  one 
state  to  another  state  is  denoted  by  the  flag  on  the  arrow  joining  the  two 
states.  Each  node  defines  the  discrete  density  probability  function  of 
finding  exactly  K cues  in  n^  glimpses. 

To  determine  the  probability  associated  with  any  specific  state, 

Pc(ng,K)  one  must  determine  all  possible  paths  from  state  (0,0)  to  the  state 
in  question.  The  probability  of  traveling  down  a specific  path  is  the  product 
of  all  the  flags  along  that  path.  The  summation  of  the  probabilities 
associated  with  all  possible  pathways  defines  Pc(ng,K).  The  probability 
of  being  in  state  (0,0)  is 

P(0,0)  - 1 

The  calculations  required  to  generate  the  density  distribution  can  be 
simplified  if  one  makes  the  observation  that  each  state  Pc(ng,K)  is  a 
function  of  at  most  two  previous  states.  For  example: 

P (3,2)  - P (2,2)  • ?(K  - 2)  + P (2,1)  * p(K  - 1)  (13) 

c c 

or  alternately 

Pc(l,0)  - Pc(0,0)  • ?(K)  (14) 
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Figure  7 Acquisition  probability  density  function 


The  expected  number  of  possible  cues  found  (E[K(ng)]  is  defined  by 


E[K(ng)] 


CO 

*/ 


K • P(ng)  ‘ dK 


(15) 


The  discrete  nature  of  this  function  allows  the  process  of  integration  to 
be  performed  by  summation.  Thus, 

n 

g 

E[K(ng)]  " 1 ' P(ng»i)  (16) 

i-1 

I 

If  all  possible  cues  have  equal  probability  of  being  the  true  cue,  the 
probability  of  finding  the  right  cue  P£  is  simply 


E[K(ng)] 

NT 


(17) 
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The  number  of  glimpses  the  user  Is  able  to  make  Is  a function  of  the 
average  glimpse  rate,  (1/tg),  the  sensor  slew  period  (tg^)  and  the  number 
of  cues  within  the  display  (ntd).  The  targets  found  within  a particular 
display  are  calculated  in  terms  of  cue  density.  Cue  density  in  an  input 
to  the  program  and  has  units  cues  per  nautical  mile  squared.  Therefore, 

the  maximum  number  of  glimpses  N the  operator  has  time  to  make  is: 

S 


N - 

g 


t + (-Si) 

8 ntd 


(18) 


The  term  t , takes  into  account  that  n. , targets  can  be  examined  without 
si  ta 

ntd 

requiring  the  user  to  slew  the  sensor. 


The  final  probability  (Pg)  to  be  calculated  is  the  probability  that  the 
target  is  within  the  sensor  field  of  view.  If  the  cues  are  uniformly 
distributed  throughout  the  search  area,  Pg  can  be  calculated  as: 

_ Area  covered  by  the  useable  footprint 
s Total  search  area 


The  area  covered  in  an  efficient  systematic  search  (Ac)  is  proportional  to 
the  total  number  of  slews  (N^)  and  the  area  of  the  footprint  determined 
useable  for  search  (A^). 


A 


c 


(20) 


The  total  search  area  (Aj)  is  assumed  to  be  a rectangle  of  length  Xj 
along  the  flight  path  and  width  2.YT.  The  total  search  area  is  defined: 


AT  - 2 • YT  • XT  (21) 

The  number  of  slews  the  operator  is  able  to  make  during  the  total  search 
process  is  defined  in  terms  of  the  platform  velocity,  search  area  length,  and 
the  amount  of  time  spent  during  each  sensor  slew.  Thus, 


T 

' IT.  - 

si 


V (t  . + n 
xv  si  g 


(22) 


Substituting  Equations  (20),  (21)  and  (22)  into  (19)  results  in 

A* 

P m i ■ ■— 

* _ ft 


Jv 
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(23) 
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Equation  (23)  shows  that  under  the  assumption  of  systematic  search  with 
a constant  depression  angle,  the  length  of  the  search  area  X-p  has  no  effect 
on  the  optimum  search  results.  The  search  area  width  on  the  otherhand  has 
a direct  effect  on  search  performance. 

In  a specific  mission  with  constant  n glimpses  per  slew,  P8  will 
decrease  as  Yp  increases.  This  decrease  in  Ps  may  be  offset  by  decreasing 
n~.  This  also  has  the  effect  of  decreasing  Pc,  but  has  an  overall  positive 
effect.  For  light  user  load  situations,  for  example  low  aircraft  velocity 
and  a "fast"  slew  time,  the  upper  bound  of  YT  may  be  the  geometric  cross- 
range capability  of  the  sensor  rather  than  the  observers  bandwidth. 

c.  VARVS  Search  Model  Evolution 

The  VARVS  search  model  has  gone  through  a series  of  revisions  in  order 
to  assure  compliance  with  experimental  data  and  to  comply  with  our  under- 
standing of  the  search  process.  Special  attention  was  taken  in  order  to 
assure  that  arbitrary  fudge  factors  were  not  introduced  into  the  model 
during  these  revisions.  The  major  revisions  to  the  model  will  be  discussed 
in  the  order  they  were  Implemented.  In  each  case  the  problem  encountered 
will  be  explained  in  terms  of  the  model  and  a brief  description  of  or 
reference  to  the  final  model  will  be  made. 


The  problems  and  changes  are: 


1)  The  original  model  was  unable  to  take  into  account  multiple  cues 
and  the  non-replacement  nature  of  cue  detection.  The  model  only  used  the 
binomial  distribution  to  calculate  the  probability  that  at  least  one  cue 
would  be  found  in  n_  glimpses. 
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This  problem  was  corrected  by  implementing  the  target  acquisition  function 
described  by  Figure  7 which  takes  into  account  both  multiple  glimpses 
and  the  effects  of  more  than  one  cue  on  the  display. 

2)  In  the  original  model  the  maximum  number  (N  ) of  glimpses  allowed 
per  sensor  position  was  determined  based  on  the  premise  that,  if  at  all 
possible,  the  total  search  area  would  be  covered.  The  lower  limit  of  n 

in  this  case  was  one.  While  this  assumption  is  sometimes  true,  it  is  not 
always  the  case.  To  assure  optimum  use  of  time,  the  process  described  by 
Equation  (2)  was  implemented. 

3)  Upon  comparing  experimental  data  with  calculated  results,  it  was 
found  that  the  calculated  prediction  for  high  cue  densities  were  sub- 
stantially lower  than  experimental  results.  Additional  comprehension  of 
the  acquisition  process  Indicated  that  the  observer  was  satisfied  to  use 
the  lower  resolution  of  the  off-axis  portion  of  the  display  in  order  to 
save  time  and  at  the  same  time  apparently  Increase  his  performance. 
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To  take  this  into  account  1)  an  average  value  of  resolution  across 
the  target  was  calculated  rather  than  the  use  of  the  on-axis  value,  2)  the 
actual  probability  contour  used  to  describe  the  useable  footprint  was 
determined  by  an  optimization  routine  and  3)  the  size  of  the  useable  foot- 
print was  used  to  determine  ntj  used  in  Equation  (18). 

The  VARVS  search  model  then  under  went  a series  of  implementation 
streamlining  to  cut  costs  and  allow  a finer  portion  of  the  footprint  (AA) 
description.  The  added  data  points  calculated  helped  give  a better  estim- 
ation of  the  effects  of  lens  distortion. 

d.  Flow  Diagram 

A simple  flow  chart  showing  the  implementation  of  the  model  is  shown 
in  Figure  8,  and  will  be  briefly  reviewed. 

First  the  ground  footprint  is  calculated  using  an  adaptation  of  the 
NVL  model.  The  footprint  is  evaluated  and  the  portion  of  the  footprint  to 
be  used  for  search  is  chosen.  This  augmented  footprint  is  characterized  in 
terms  of  length,  width,  and  average  number  of  line  pairs  across  the  target. 
These  characteristics  are  used  to  calculate  the  time  available  for  search 
(ta),  and  the  probability  of  detecting  the  cue  in  a single  glimpse  given  there 
is  a single  cue  and  the  observer  is  looking  at  the  cue  (P'). 

The  first  optimization  to  be  made  is  to  allocate  the  time  available 
between  pure  display  search  and  slew  time.  If  the  user  spends  too  much 
time  looking  at  one  portion  of  the  total  search  area  he  is  apt  to  leave 
large  portions  of  the  search  area  unsearched.  Alternatively,  if  the  observer 
wastes  too  much  time  slewing  the  sensor,  he  cuts  down  on  the  amount  of  time 
he  has  for  search.  The  indicated  target  acquisition  function  (TAF)  loop 
determines  the  optimum  combination  of  glimpses  per  slew  and  area  coverage 
as  it  pertains  to  each  situation. 

When  an  optimum  solution  is  found,  the  probabilities  of  target  detection 
(Pd)  and  recognition  (Pr)  are  calculated.  The  probability  of  mission  success 

is  evaluated  in  terms  of  the  probability  of  detecting  the  target  (pDET) 

PDET  ' Pc  • Pd  ‘ P.  (25> 

and  the  probability  of  recognizing  the  target  (P?) 


p ■ p • P 
rT  DET  r 

The  probabilities  of  mission  success  are  evaluated  and  a new  strategy  of 
search  may  be  chosen  if  the  possibility  of  better  results  is  indicated. 

A new  strategy  may  either  Increase  or  decrease  the  size  of  the  useable 
footprint.  For  example,  decreasing  the  useable  footprint  size  will  increase 
PT  and  force  the  product  of  Pc  and  P8  to  decrease.  When  an  adjustment  is 
made  in  the  size  of  the  footprint,  the  time  allocation  optimization  must  be 
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reevaluated  and  then  PpET  and  Pr  are  recalculated.  The  loop  continues  until 
the  model  operator  determines  that  the  optimum  has  been  found.  This  shows 
the  adaptiveness  that  must  be  present  in  the  model  and  in  real  search  to 
arrive  at  optimum  results. 

The  TAF  loop  is  implemented  on  a digital  computer  while  the  decision 
to  try  a different  portion  of  the  footprint  is  left  to  the  program  operator. 
The  major  cost  in  computer  time  is  spent  in  the  generation  of  the  ground 
footprint.  Because  of  the  lens  distortion  function  a new  footprint  must 
be  generated  at  each  depression  angle,  altitude  or  velocity.  This  cost 
in  computer  time  to  generate  new  footprints  becomes  the  limiting  factor  as 
far  as  cost  is  concerned. 


om-min 


Figure  9 Moving  terrain  teat  apparatus 


3.  EXPERIMENT  DEVELOPMENT 

The  experiment  consisted  of  a 10  ft.  x 20  ft.  treadmill  which  gave  a 
10  ft.  by  40  ft.  simulated  terrain  map,  a camera,  and  remote  display 
stations.  The  treadmill  and  camera  are  shown  in  Figures  9 and  10.  A 
600/1  scale  was  selected.  This  provided  useable  coverage  area  of  6000  ft. 
by  24,000  ft.  (4  sq.  miles).  At  a simulated  velocity  of  1,000  ft/sec.  the 
mission  time  is  24  seconds.  The  camera's  optical  axis  was  located  3.5  ft. 
above  the  terrain  map  which  represents  a scale  altitude  of  2,100  ft.  Thus 
the  sensor  could  operate  over  scale  slant  ranges  varying  from  12,800  ft.  to 
about  2,200  ft.  when  the  sensor  is  located  one  foot  from  the  treadmill. 
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Two  areas  of  difficulty  were  encountered  with  the  experiment.  The 
first  was  in  getting  the  terrain  map  of  the  treadmill  to  run  straight  and 
true.  This  was  corrected  by  installing  a center  spreader  between  the 
rollers  of  the  treadmill.  This  caused  the  rollers  to  bow  outwards  at 
their  center  while  maintaining  the  constant  linear  surface  velocity  that 
is  necessary  to  prevent  wrinkles  from  forming  in  the  map  (this  is  also  the 
reason  the  conventional  solution  to  belt  tracking  problems,  domed  rollers, 
could  not  be  used).  The  other  problem  area  was  in  selecting  the  best  type 
of  cues  and  targets  for  the  experiment.  Since  the  analytic  model  is  based 
on  the  observer  visually  examining  a number  of  cues  to  determine  which  one 
contains  the  target,  we  decided  the  cues  could  be  identical.  Since  many 
cues  wgre  required,  this  greatly  reduced  the  effort  required  to  fabricate 
the  map.  We  found  it  very  difficult,  however,  to  find  a cue  which  did  not 
reveal  the  target  because  of  contrast  variation.  For  example,  the  target 
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Figure  1 0 Camera  tracking  targets  on  terrain 
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must  have  a critical  dimension  of  1/5  the  cue  dimension.  If  one  uses  a 
square  or  circle  as  a cue  with  a smaller  contrasting  square  or  circle  as 
a target,  the  fact  that  the  cue  contains  a target  is  instantly  visible 
because  this  cue  now  has  a noticeably  lower  contrast  than  the  other  cues. 

This  led  us  to  use  Landholt  C ring  as  cues  and  targets.  The  Landholt  C 
ring  shown  in  Figure  11,  is  well  accepted  for  acquisition  studies.  In 
our  experiment  solid  rings  were  used  except  for  cues  with  a gap  which 
represent  the  target.  The  gap  occupies  a very  small  portion  of  the  ring 
and  therefore  produces  no  noticeable  contrast  loss. 

Implementing  the  C ring  for  our  experiment  proved  much  more  difficult 
than  anticipated.  A circle  on  the  terrain  map  appears  as  an  ellipse  on 
the  viewing  surface.  This  makes  the  gap  have  variable  angular  subtense 
depending  on  its  location  on  the  C ring  and  the  viewing  geometry.  We 
corrected  this  problem  by  locating  the  gap  on  a side  of  the  ring  and  by 
making  the  gap  three  dimensional.  This  created  another  problem,  however. 

Due  to  shading  of  the  ring  by  the  3D  gap,  a noticeable  contrast  variation 
in  the  target  ring  relative  to  the  closed  rings  were  created.  To  correct 
this,  the  C rings  were  made  elliptical  and  oriented  so  that  they  appear 
circular  and  have  uniform  crossection  on  the  viewing  screen.  Since  the 
viewing  geometry  continually  varies  over  the  terrain  map  we  chose  to 
optimize  geometry  for  what  we  thought  to  be  the  nominal  acquisition  dlstance- 
between  6,000  and  8,000  ft.  downrange.  This  ring  configuration  is  shown 
in  Figure  12.  We  determined  experimentally  that  this  ring  worked  quite 
well  over  the  entire  12,000  x 6,000  scale  feet  viewing  area.  At  short  ranges 
the  elliptical  shape  was  apparent  but  the  rings  occupied  so  large  an  angular 
subtense  that  it  caused  no  distraction. 

Since  the  threshold  resolution  of  the  VARVS  is  in  the  area  of  one 
mllliradlan,  we  decided  to  scale  the  C rings  to  provide  a one  milliradian 
target  subtense  at  the  maximum  slant  range.  This  would  make  target  detection 
possible  anywhere  on  the  map  and  make  search  dynamics  the  dominant  factor. 

For  convenience  in  fabrication,  a 0.25  inch  target  dimension  was  chosen. 

This  yields  at>  angular  subtense  of  0.99  mr  at  the  maximum  slant  range.  The 
other  ring  dimensions  were  scaled  from  Figure  12. 

In  the  actual  experiment  the  solid  cues  (rings)  were  placed  randomly 
on  the  map  with  the  required  average  cue  density  as  shown  in  Figure  13.  The 
cue  containing  the  target  was  placed  on  the  map  randomly  at  the  start  of 
each  mission.  The  observer's  task  was  therefore  to  find  the  target  before 
it  was  overflown.  The  exact  procedures  will  be  discussed  later. 

Another  detail  that  had  to  be  worked  out  before  the  experiment  could 
proceed  was  the  number  of  runs  or  missions  required  to  establish  the 
statistical  significance  of  the  experiment  data.  An  analysis  was  run  using 
the  Chebyshev  Inequality.  This  analysis,  presented  in  Appendix  C yielded 
the  equation: 

»>  (27) 

Y • e 
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Figure  1 3 Cue  density  scaling 


n is  the  number  of  data  points  required 
p Is  the  average  probability  of  detection  expected 
e is  the  desired  tolerance  interval 

Y is  the  probability  data  points  will  fall  out  of  the  tolerance 
band  +e 


This  equation  is  plotted  on  Figure  14  for  10%  and  5%  tolerance  levels. 

For  the  target  scenario  chosen,  the  preliminary  performance  analysis  indicated 
about  80%  probability  of  target  detection  should  be  expected.  Figure  14 
shows  about  100  data  points  would  be  required  to  establish  Pj  within  10% 
while  400  points  would  be  required  to  establish  within  5%.  Considering 
the  great  number  of  conditions  that  must  be  run,  we  elected  to  use  the  10% 
tolerance  level. 

Since  many  data  points  must  be  collected  for  each  set  of  input  para- 
meters, an  accurate  estimate  of  the  parameter  ranges  to  be  studied  was 
required.  As  a minimum  we  felt  we  should  run  four  depression  angles, 
three  speeds,  two  cue  densities,  and  two  offsets  for  each  of  the  four 
sensor  types.  In  addition  we  wanted  to  run  the  VARVS  with  both  hand  and 
head  control  using  the  ONR  dome  projection  system  and  hand  control  with  the 
USAF  direct  view  display.  Altitude  was  not  Included  because  of  the 
difficulty  of  Implementing  changes  to  the  test  set-up  and  these  effects 
would  be  purely  geometrical  and  are  easily  predicted.  The  actual  test  plan 
is  shown  in  Figure  15. 

The  actual  test  was  run  as  follows.  A test  conductor  was  stationed 
at  the  lower  end  of  the  treadmill  adjacent  to  the  camera  as  shown  in 
Figure  16.  The  test  conductor  had  a two  way  communication  link  with  the 
observer  who  was  at  the  display  station  in  another  room.  A third  man,  the 
map  controller  was  located  under  the  treadmill  as  shown  in  Figure  16.  The 
map  controller's  job  is  to  continually  move  the  cues  and  targets  on  the  map 
so  the  observer  does  not  memorize  the  pattern.  The  test  director  tells  the 
observer  what  depression  angle  to  use  and  checks  visually  to  see  that  the 
camera  is  properly  oriented  downrange.  He  then  gives  a start  command  to  the 
map  controller.  After  a variable  time  delay  (0  to  15  sec.)  the  controller 
places  a cue  on  the  map  which  contains  the  target.  Meanwhile  the  observer 
searches  the  map  for  the  target.  If  the  observer  finds  the  target,  he 
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Figure  1 6 Test  plan 
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Figure  1 6 Experimental  test  setup 


notifies  the  test  conductor  who  checks  the  camera  orientation  and  auxiliary 
monitor  to  see  that  it  is  truly  pointing  at  the  target.  If  it  is  a valid 
acquisition  the  test  is  terminated.  Otherwise  the  test  continues  until 
the  conductor  sees  the  target  has  been  overflown.  Generally  ten  missions 
were  run  before  the  observer  is  given  a rest.  Observers  were  changed  after 
four  sets  of  ten  missions  each.  After  all  missions  were  completed  with 
one  type  of  sensor,  the  sensor  was  converted  to  another  type.  In  order 
to  keep  all  system  parameters  constant,  only  the  sensor  lens  was  changed. 
For  the  VARVS,  the  non-linear  lens  was  used.  For  linear  sensor  tests,  the 
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non-linear  lens  was  removed  from  the  VARVS  sensor  and  replaced  with  a 
35-105  mm  F/3.5  Vivitar  zoom  lens.  This  lens  was  adjusted  for  the 
required  f ield-of-view  and  the  next  set  of  data  was  run. 

4 . MODEL  VALIDATION 

Two  efforts  on  model  validation  were  accomplished.  The  first  is 
a theoretical  comparison  of  results  predicted  by  our  new  linear  sensor 
model  with  those  predicted  by  experimentally  validated  existing  models. 

The  other  is  the  comparison  of  both  VARVS  and  linear  sensor  theoretical 
and  experimental  results. 

To  compare  with  existing  models,  the  dynamic  performance  of  a 20°  TV 
was  computed  by  both  Bailey's  model  (Reference  6)  and  by  our  new  linear 
sensor  model.  The  results  are  shown  on  Figures  17,  18,  and  19  for 
various  depression  angles.  At  lateral  coverages  greater  than  4000  feet, 
results  are  in  good  agreement  if  a clutter  factor  of  one  is  used  in  the 
Bailey  model  and  a cue  density  of  about  60  cues/nm^  is  used  in  our  new 
model.  At  smaller  offsets  the  new  model  predicts  better  performance  than 
Bailey's  model.  This  is  due  entirely  to  peripheral  cueing  which  becomes  v 
more  effective  when  more  search  time  is  available  at  small  lateral  coverages. 
Considering  the  empirical  nature  of  Bailey's  clutter  factor  and  difficulty 
in  relating  our  cue  density  to  the  real  world,  we  believe  the  agreement 
between  models  is  fairly  good.  The  main  observation  to  be  made  here  is 
that  since  Bailey's  clutter  factor  of  one  represents  a fairly  simple  scene 
such  as  a desert  and  this  corresponds  to  a cue  density  of  60,  a rather  high 
cue  density  in  the  order  of  200-300  cues/nm^  is  required  to  represent 
real  world  scenes. 

The  experimental  results  are  next  compared  to  the  theoretical  predic- 
tions. Input  parameters  for  the  theoretical  predictions  are  shown  on 
Figure  20.  The  math  models  used  were  the  static  and  dynamic  models  described 
in  Sections  3.1  and  3.2.  The  comparison  is  shown  in  Figures  21  through  29. 
The  input  parameters  were  based  on  system  performance  measurements  i.e,, 
theoretical  parameters  were  adjusted  to  correlate  with  the  experimental 
numbers.  These  performance  numbers  were  checked  at  the  beginning  and  end 
of  each  test  day  using  the  test  target  at  the  far  end  of  the  treadmill 
(See  Figure  9 or  10). 

One  factor  that  was  observed  early  in  the  experimental  effort  was  that  a 
large  difference  exists  between  resolving  C ring  targets  and  the  3 bar  test 
targets.  In  order  to  study  this,  a C ring  test  target  was  placed  next  to 
the  3 bar  target  - See  Figure  9 or  10.  The  system  was  defocused  until  a C 
ring  was  just  resolved.  Then  the  frequency  of  the  3 bar  pattern  that  was 
also  just  resolved  was  noted.  Regardless  of  the  C ring  size,  we  found  that 
the  3 bar  pattern  having  one  half  bar  width  equal  to  the  C ring  gap  could 
be  resolved  also.  Since  the  model  was  developed  and  adjusted  to  predict 
3 bar  resolution,  we  had  to  use  a target  scale  size  of  25  ft.  to  predict 
experimental  results  which  used  a C ring  gap  of  scale  12.5  ft.  Figures 
21  through  24  show  comparisons  for  the  four  system  types,  the  VARVS  and 
linear  sensor  with  three  f ields-of-view.  On  each  of  these  curves  the 
theoretical  prediction  from  the  new  model  are  plotted  as  solid  lines.  The 
experimental  points  are  shown  as  averages  of  the  100  data  points  and  the 
associated  standard  deviation  in  the  data.  In  general  one  will  observe 
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Probability  of  Detecting  a True  Cue  Probability  of  Detecting  a True  Cue 


Lateral  Coverage  ■ 1000  ft  onmiMi 


Figure  17  Probability  of  detection  of  TV  with  20°  FOV 
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Figure  1 9 Probability  of  detection  of  TV  with  20°  FOV 


Probability  of  Target  Detection  Probability  of  Target  Detection 


Downrange  Pointing  Location  of  Optical  Axis  - 1000  ft 


Figure  24  Comparison  of  35.5°  FOV  TV  experimental  to  theotetical  data 
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• Curves  represent  theoretical  data 

• Points  represent  experimental  data 
(average  and  one  sigma  excursion) 
plotted  at  the  average  acquisition  range. 
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Figure  25  Effect  of  cue  density  on  the  perfomance  of  the  RVS 
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that  the  upper  standard  deviation  is  very  close  to  the  theoretical  data. 
This  is  to  be  expected  because  the  model  assumes  consistent  and  mflvtfmtm 
human  performance  at  all  times.  It  was  obvious  throughout  the  test  that 
this  was  not  the  case  with  any  observer.  His  performance  would  vary  with 
such  things  as  time  of  day,  test  duration,  etc.  It  seems  logical  therefore 
that  the  best  of  the  experimental  data  should  match  the  theoretical 
predictions.  Some  thought  was  given  at  this  point  in  the  study  to 
including  observer  variability  in  the  model.  This  was  not  done  because, 

(1)  the  comparative  nature  of  the  parametric  study  should  make  results 
independent  of  this  factor  and  (2)  the  model  had  undergone  numerous 
iterations  to  achieve  this  level  of  agreement  and  time  was  running  out 
to  make  further  changes.  Model  validation  for  cue  density  and  velocity  are 
given  in  Figures  25  through  28. 

A comparison  of  the  h«.nd  controlled  mode  to  the  head  controlled 
mode  is  shown  on  Figure  29.  The  reason  for  lower  performance  for  the 
head  control  mode  is  the  additional  servo  lag  caused  by  the  projector 
gimbals.  The  system  slew  time  increases  from  0.5  sec  to  0.7  sec.  when 
the  projector  servo  is  activated. 

In  general  velocity  and  cue  density  trends  are  correct.  However,  as 
noted,  insufficient  data  points  were  taken  in  many  cases.  It  is  disturbing 
that  the  high  cue  density  data  for  the  VARVS  exceeds  the  theoretical  data. 

We  wanted  to  run  a 40  cue/nm^  experiment  to  see  if  the  trend  continues  but 
time  did  not  permit  this.  No  data  for  lateral  coverage  is  presented  because 
for  the  range  of  coverages  possible  with  the  test  setup  (0-1500  feet) 
experimental  and  theoretical  probability  data  were  100%. 

The  experiment  using  the  direct  view  display  was  terminated  after 
several  missions  because  all  observers  complained  that  it  was  very 
fatiguing  to  use. 
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SECTION  IV 


PARAMETRIC  ANALYSIS 


In  this  section  the  performance  of  the  IR  VARVS  is  compared  to  a con- 
ventional IR  TV  system  using  the  analysis  models  developed  and  identified 
earlier.  Parameters  of  cue  density,  altitude,  velocity  and  target  offset 
from  the  flight  path  are  pertubated  to  measure  the  effectiveness  of  the  two 
systems.  This  effectiveness  then  can  be  related  to  the  particular  mission 
type  best  suited  for  the  sensor.  It  was  determined  that  in  addition  to  the 
above  parameters,  the  IR  VARVS  was  sensitive  to  slew  time  so  this  was  also 
retained  in  the  parametric  analysis. 

The  infrared  target  acquisition  models  consists  of  two  static  sensor 
performance  models  and  a dynamic  model.  The  static  models  are  the  NVL  model 
(Reference  4)  for  the  linear  sensor  and  the  MCAIR  IR  VARVS  model  which  is  an 
adaptation  of  the  NVL  model.  This  model  was  developed  in  prior  efforts 
(Reference  2 and  3)  and  is  included  as  Appendix  D.  System  parameters  used 
in  modeling  the  IR  VARVS  and  three  conventional  FLIRS  are  shown  in  Figure  30. 
The  FLIR  scanners  are  second  generation  types  with  488  resolution  elements, 
parallel  array  with  a nine  element  Time  Delay  Integration  (TDI)  which  was 
shown  in  the  Reference  2 report  to  be  optimum  for  the  IR  VARVS.  The  optical 
quality  for  the  non-linear  lens  was  also  taken  from  that  study.  Display 
quality  for  the  IR  VARVS  in  the  model  is  based  on  state  of  the  art  in  optical 
fabrication  and  is  somewhat  better  than  the  experimental  equipment  (50  p blur 
compared  to  105  p) . FLIR  system  parameters  are  believed  to  be  representative 
of  second  generation  state  of  the  art  systems. 

The  approach  and  data  reduction  techniques  used  in  parametric  analysis 
were  described  in  Section  II.  A typical  sequence  in  data  reduction  for  the 
reconnaissance  mission  is  presented  below.  The  IR  VARVS  and  20°  FLIR  cases 
for  the  nominal  parameters  will  be  used  for  this  purpose. 


Parameter 

IR  VARVS 

10°  FLIR 

20°  FLIR 

40°  FLIR 

Number  of  Scan  Lines 

525 

525 

525 

525 

Light  Level 

<fL) 

10 

10 

10 

10 

Lens  Blur 

(microns) 

16 

16 

16 

16 

Display  Blur 

(microns) 

50 

50 

50 

50 

Detector  Site 

(in.) 

0.001 

0.001 

0.001 

0.001 

Bandwidth 

(Me) 

- 

- 

- 

- 

F/No. 

2.0 

2.0 

2.0 

2.0 

Altitude  (Nominal) 

(ft) 

1500 

1500 

1500 

1500 

Velocity  (Nominal) 

(ft/sec) 

840 

840 

840 

840 

Slew  Period  (Nominal) 

(sec-1) 

0.3 

0.3 

0.3 

0.3 

Glimpse  Time 

(sec) 

0.25 

0.25 

0.25 

0.25 

Resolution 

(mr/cycle) 

0.925 

0.436 

0.877 

1.81 

Focal  Length 

(mm) 

50.8 

101.6 

50.4 

74.4 

opra-ina-i 

Figure  30  System  peremeters 


On  Figure  31  the  model  probability  of  detection  outputs  are  plotted  as 
a function  of  the  lateral  coverage  for  various  cue  densities.  Lateral 
coverage  values  for  90%  probability  of  detection  are  extracted.  These  plus 
those  from  similar  plots  for  other  depression  angles  are  used  to  plot  Figure  32. 
Using  this  curve,  optimum  depression  angle  can  be  determined  based  on  lateral 
coverage.  In  the  sample  case  the  optimum  depression  angle  is  about  14°.  The 
lateral  coverage  at  this  depression  angle  will  be  used  as  a figure  of  merit  of 
the  system.  For  the  conventional  FLIR  a similar  procedure  is  used  except  the 
computer  selects  the  90%  probability  automatically.  Therefore,  the  first  step 
is  to  make  a depression  angle  plot  shown  in  Figure  33.  Here  some  explanation 
of  the  shape  of  this  plot  is  in  order.  First  one  must  remember  that  the 
depression  angle  is  depression  of  the  sensor's  optical  axis  and  bears  no 
strict  relationship  to  the  actual  range  where  the  observer  elects  to  search. 

As  discussed  in  Appendix  C,  the  observer  will  use  only  the  portion  of  the 
display  that  presents  suitable  quality  to  insure  a reasonable  probability 
of  target  acquisition.  It  is  therefore  possible  that  the  same  "useable" 
footprint  will  be  available  to  the  observer  over  a range  of  depression 
angles.  When  the  depression  angle  outside  of  the  optimum  range  is  chosen, 
the  useable  search  area  will  be  cut  off  by  the  sensor  format  edge.  This 
explains  the  sudden  rise  and  drop  in  performance  at  both  low  and  high  depres- 
sion angles.  The  performance  curves  indicate  a subtle  increase  in  acquisition 
performance  as  the  depression  angle  goes  from  20°  to  3°.  This  trend  is 
especially  apparent  for  the  lower  target  densities.  The  increase  in 
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Figure  31  IR  VARVS  probability  of  detection  for  14°  depression  angle 
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Figure  32  IR  VARVS  dynamic  coverage  capability  vs  depression  angle 


performance  with  shallower  depression  angles  is  primarily  a function  of  the 
increased  sensor  ground  coverage.  The  larger  footprint  size  increases  the 
time  available  for  search  and  results  in  less  sensor  slewing.  The  longer 
slant  range  will  tend  to  decrease  ground  resolution  but  in  the  20°  F0V  FLIR 
case  shown,  the  net  effect  results  in  an  increase  in  the  performance  measure. 
This  phenomenon  will  be  used  later  to  select  the  minimum  weapons  delivery 
range.  The  effect  of  depression  angle  is  not  as  pronounced  as  the  field- 
of-view  gets  wider. 

As  with  the  IR  VARVS,  the  peak  lateral  coverage  is  selected  from  Figure  33 
and  used  as  the  figure  of  merit  of  the  conventional  FLIR.  The  relative  per- 
formance plots  are  then  made  by  dividing  the  VARVS  peak  lateral  coverage  by 
the  peak  FLIR  lateral  coverage.  The  results  for  the  sample  case  is  shown 
on  Figure  34.  If  the  ordinate  goes  above  one,  the  IR  VARVS  is  better 
chan  the  FLIR. 

The  weapon  delivery  plots  require  some  additional  effort.  Here  the 
figure  of  merit  is  peak  downrange  acquisition  distance  and  lateral  coverage 
product.  With  the  VARVS,  the  search  process  is  forced  to  be  near  the  optical 
axis  because  of  acuity  fall-off.  Therefore  the  acquisition  range  is  very 
near  the  range  given  by  the  depression  angle.  The  actual  weapons  delivery 
range  can  be  anywhere  within  the  probability  contour  chosen  for  search.  This 
contour  is  described  in  Section  2 and  has  nominally  been  found  to  be  the  90% 
target  detection  contour.  The  minimum  weapons  delivery  range  therefore 
corresponds  to  the  downrange  distance  to  the  aft  edge  of  the  90%  contour. 

Since  Figure  32  shows  that  lateral  coverage  is  essentially  constant  for  the 
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Figure  33  20°  FLIR  dynamic  coverage  capability  vs  depression  angle 
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Figure  34  Relative  lateral  coverage  performance  of  VARVS  end  20°  FOV  FLIR 
for  reconnaissance,  sample  case 


IR  VARVS,  the  greatest  range  and  lateral  coverage  product  will  be  at  the 
minimum  useable  depression  angle.  The  distance  calculated  from  this  angle 
and  adjusted  to  correspond  to  the  aft  edge  of  the  footprint,  multiplied  by 
the  lateral  coverage  is  the  figure  of  merit  of  the  IR  VARVS  sensor. 

For  the  conventional  FLIR,  calculation  of  acquisition  range  is  some- 
what more  difficult  because  lateral  coverage  is  a function  of  depression 
angle.  Fortunately,  the  peak  lateral  coverage  always  occurs  at  minimum 
depression  angle.  The  maximum  acquisition  range  is  the  downrange  distance 
of  the  aft  edge  of  the  useable  footprint  at  minimum  depression  angle.  This 
also  corresponds  to  maximum  lateral  coverage.  To  be  compatible  with  the  IR 
VARVS*  figure  of  merit,  the  product  of  this  range  and  the  lateral 
coverage  is  defined  as  the  figure  of  merit  of  the  FLIR.  For  the  sample 
case,  the  values  are  as  shown  in  Figure  35.  Figure  36  shows  the  curve  for  the 
sample  case. 

1.  RESULTS 

The  results  comparing  IR  VARVS  to  the  three  FLIR  configurations  are 
shown  in  Figures  37-44 . It  should  be  noted  that  all  results  of  the  para- 
metric analyses  are  at  the  target  detection  level  rather  than  recognition. 

This  is  because  recognition  was  found  to  be  either  impossible  due  to  lack 
of  resolution  or  not  practical  because  of  the  short  range  required.  All  of 
the  systems  studied  would  need  auxiliary  optics  to  achieve  a good  recognition 
capability.  The  linear  FLIR s would  require  a narrower  f ield-of-view  or  zoom 
optics  while  the  VARVS  would  need  an  auxiliary  optical  system  with  five  to 
ten  times  more  magnification.  Such  a system  is  discussed  in  Section  6. 
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Figure  36  Maximum  acquisition  range  for  weapons  delivery 
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Figure  36  Relative  performance  of  IR  VARVS  and  20°  FOV  FLIR  for 
weapons  delivery,  sample  case 


Altitude  1500  ft 
Velocity  ■ 840  ft /tec 
Cue  tire  60  ft 
Target  tire  -12  ft 
Cue  dentuy  40  cuet/nm^ 
Glimpte  time  0.25  tec 


oert  mtu 


Figure  37  Relative  lateral  coverage  performance  of  IR  VARVS  and  FLIR  for 
reconnaissance  showing  sensitivity  to  slew  time 
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Figure  40  Relative  lateral  coverage  of  IR  VARVS  and  FUR  for 
reconnaissance  showing  sensitivity  to  velocity 
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Figure  41  Relative  performance  of  IR  VARVS  and  FLIR  for  weapons  delivery 
showing  sensitivity  to  slew  time 
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Figure  42  Relative  performance  of  IR  VARVS  and  FUR  for  weapons  delivery 
showing  sensitivity  to  cue  density 
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Figure  43  Relative  performance  of  IR  VARVS  and  FUR  for  weapons  delivery 
showing  sensitivity  to  altitude 
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Figure  44  Relative  performance  of  IR  VAR  VS  and  FUR  for  weapons  delivery 
showing  sensitivity  to  velocity 


Figures  37-44  show  performance  sensitivities  to  the  IR  VARVS  system  parameter 
of  slew  time  and  the  mission  parameters  of  cue  density  (clutter),  altitude, 
and  velocity.  T^e  nominal  parameters  used  in  these  plots  are  0.3  sec.  slew 
time,  40  cues/nm  cue  density,  1,500  ft.  altitude  and  840  ft. /sec.  velocity. 
These  were  selected  because  they  were  though  to  be  realistic  except  for  cue 
density.  A higher  value  would  have  been  more  better.  However,  computational 
problems  were  encountered  which  would  have  prevented  plotting  altitude  and 
velocity  sensitivity  plots.  The  cue  density  sensitivity  plots  show  that  the 
IR  VARVS  would  perform  better  compared  to  the  FLIR  at  higher  cue  density 
values  than  at  the  40  cues/nm2  value,  but  the  effect  is  rather  small  (Figures 
34,  36,  38  and  42). 

Looking  first  at  the  reconnaissance  plots,  Figure  37  shows  the  effect 
of  varying  the  VARVS  slew  time.  The  VARVS  sensitivity  to  slew  time  is  very 
great.  However,  it  does  not  approach  the  performance  of  the  FLIR  until  slew 
time  approaches  zero.  The  VARVS  simply  uses  up  too  much  of  the  available 
time  to  slew  from  cue  to  cue.  This  excess  slewing  is  required  to  move  cues 
into  the  higher  resolution  portion  of  the  display.  The  linear  FLIR  is  more 
time  efficient  because  slewing  is  only  required  to  achieve  the  search  area 
coverage  i.e.,  tar  can  be  equally  detected  anywhere  within  its  field- 
of-view.  The  slev  used  in  the  experiment  was  about  0.5  seconds  but 

the  time  response  t -ulpment  could  be  Improved  to  decrease  the  slew  time 
to  0.3  seconds,  this  value  was  used  as  nominal  for  the  parametric  analyses. 
Theoretically,  slew  time  could  be  further  reduced  to  0.1  or  0.05 
seconds.  However,  eye  control  rather  than  head  control  would  be  required 
to  realize  the  benefits  of  this  Improvement. 


44 


From  the  data  of  Figure  37,  one  can  conclude  that  the  basic  head  con- 
trolled IR  VARVS  is  limited  to  50-80%  of  conventional  linear  FLIR  perfor- 
mance. With  eye  control  this  could  be  improved  to  the  80  to  100%  range. 

The  remainder  of  the  reconnaissance  plots,  Figures  38,  39  and  40  show 
sensitivities  of  the  comparison  to  mission  parameters.  The  cue  density 
data,  Figure  38,  shows  that  {he  IR  VARVS  performance  improves  relative  to 
the  FLIR  as  cue  density  or  clutter  is  increased.  This  trend  occurs  because 
at  higher  cue  densities,  multiple  cues  appear  within  the  portion  of  the  IR 
VARVS  footprint  where  target  detection  is  likely.  This  means  that  more 
cues  can  be  inspected  for  the  target  without  slewing  the  sensor.  The 
altitude  and  velocity  plots,  Figure  39  and  40,  respectively,  show  that  the 
nominal  parameters  of  1,500  ft.  and  840  ft/sec.  are  very  near  the  optimum 
and  that  the  IR  VARVS  performance  reduction  on  either  side  of  these  points 
are  small  but  significant. 

The  weapon  delivery  plots.  Figures  41,  42,  43  and  44,  show  how  the 
product  of  downrange  acquisition  distance  and  lateral  coverage  vary  with 
the  same  sensor  and  mission  parameters  that  were  used  on  the  reconnaissance 
curves.  Here  the  results  are  somewhat  different.  The  performance  for 
weapon's  delivery  are  more  dependent  on  sensor  fleld-of-view  as  shown 
because  the  curves  for  the  10°,  20°,  and  40°  systems  are  more  widely  spaced. 
This  is  a direct  result  of  the  FLIR  resolution  capability  which  affects  the 
probable  target  acquisition  range.  In  the  weapon  delivery  case,  the  IR  VARVS 
performance  exceeds  that  of  the  40°  FLIR  for  the  nominal  mission  parameters. 

On  the  other  hand,  the  VARVS  has  only  40%  of  the  performance  of  a 20°  FLIR 
which  has  similar  resolution  capability.  The  VARVS  is  even  less  effective 
than  a 10°  FLIR.  This  is  because  lateral  coverage  capabilities  are  comparable 
(Figure  37-40)  while  the  acquisition  range  is  much  greater  for  the  higher 
resolution  10°  TV. 
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SECTION  V 


OPTIMUM  PARAMETER  DISCUSSION 


The  thorough  knowledge  of  the  VARVS  gained  from  the  theoretical  and 
experimental  studies  described  herein  makes  the  path  clear  to  achieve  optimum 
performance  from  the  system.  The  rationale'  develops  as  follows. 

The  analyses  shows  that  two  routes  exist  to  Improve  VARVS  performance, 
by  reducing  slew  time  and/or  Increasing  the  size  of  the  high  resolution  area 
of  the  final  displayed  Image.  Decreasing  slew  time  will  make  search  more 
efficient  by  making  more  slews  possible  while  Increasing  the  high  acuity 
area  requires  fewer  slews.  There  appears  to  be  a tradeoff  here.  However, 
reference  to  the  parametric  analysis.  Section  3,  provides  some  Interesting 
observations.  By  reducing  the  slew  time  to  zero.  Figures  37  and  41  show 
that  VARVS  performance  equals  the  20°  linear  sensor.  When  one  realizes  both 
sensor  types  have  identical  on-axis  resolution  performance,  broadening  the 
on-axis  acuity  will  increase  performance  until  the  point  of  uniform  acuity 
across  the  entire  visual  field  Is  reached.  We  would  then  have  a conventional 
linear  sensor.  If  a linear  sensor  has  target  acquisition  performance  equal 
to  the  VARVS,  then  the  VARVS  advantage  must  lie  ei'sewhere,  per naps  for  flight 
control  and  navigation  where  the  advantages  of  a wide  f leld-of-vlew  is  clear. 

If  one  accepts  the  value  of  the  wide  field,  the  alternate  route  of 
improving  target  acquisition  performance  of  the  VARVS  must  be  persued.  The 
data  shows  that  slew  time  must  be  reduced  to  a fraction  of  the  glimpse  time 
in  order  to  achieve  performance  comparable  to  a conventional  linear  sensor. 

Slew  times  on  the  order  of  0.01  to  0.03  sec.  would  be  required  by  eye  control 
to  achieve  the  required  performance.  This  would  be  a major  technical  task. 
Perhaps  a better  way,  made  possible  by  advancements  in  solid  state  array 
technology,  would  be  to  use  a wide  field  linear  lens  with  a focal  plane  array 
capable  of  supporting  high  resolution  anywhere  within  its  f ield-of-view.  The 
variable  acuity  would  then  be  obtained  electronically  by  sampling  the  array 
in  the  proper  manner.  The  high  speed  eye  controlled  slewing  could  then  be 
handled  electronically  in  a field  time  of  0.017  sec.  With  this  capability, 
the  mechanical  sensor  slewing  would  then  only  have  to  follow  much  slower  head 
motion.  The  focal  plane  array  required  by  this  approach  would  be  about 
10000  x 10000  elements.  Detector  elements  would  have  to  be  very  small,  no 
more  than  0.0003  inches  to  produce  a reasonable  array  size  of  3 Inches  square. 
The  display  would  require  similar  performance  (10,000  lines)  and  be  helmet 
mounted  to  support  this  approach.  Such  a system  must  await  sensing  and  emitting 
array  state  of  the  art  advancements. 


46 


SECTION  VI 


DUAL  FIELD  FEASIBILITY  STUDY 


Analyses  of  Section  V have  documented  IR  VARVS  performance  to  be  In 
the  0. 8-1.0  mllliradlan  area.  This  Is  suitable  for  target  detection  but 
Is  unuseable  for  target  recognition  tasks.  For  the  recognition  task, 
according  to  the  Johnson  Criteria,  about  Sx  magnification  would  be  required 
so  that  the  target  subtends  an  adequate  number  of  lines. 

Since  the  wide  angle  optical  quality  is  not  sufficient  to  support  this 
performance,  a separate  optical  system  must  be  considered.  Such  an  approach 
that  uses  a common  aperture  Is  shown  in  Figure  A 5.  Here  the  front  element 
of  the  wide  angle  lens  is  comnon  with  narrow  f ield-of-view  high  resolution 
optical  system.  The  wide  field  optical  system  is  that  described  in 
Reference  2 which  Is  believed  to  be  the  optimum  design  for  an  IR  VARVS. 

The  low  optical  power  of  the  large  front  element  to  nearly  paraxial  rays 
makes  possible  Its  use  as  the  first  element  of  a long  focal  length  high 
resolution  optical  system,  providing  that  the  small  optical  deviations 
introduced  by  the  aspheric  element  can  be  corrected  outside  of  the  optical 
path  occupied  by  the  wide  angle  optical  system.  This  is  the  function  of 
the  corrector  element  shown  on  Figure  45.  The  feasibility  of  this  corrector 
and  remainder  of  the  narrow  field  optical  system  is  described  below. 


- Axial  Edge  Ray  for  Wide  Angle  Mode 


Selector  Mirror 


■ Axial,  0.45  of  Ref  Pupil  Height,  Narrow 
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- Chief  Ray  of  Edge  Field 
Wide  Angle  Mode 
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Figure  46  Common  apsrture  VARVS 


1 . INPUT  PARAMETERS 


The  first  step  to  establish  the  feasibility  of  a narrow  field  optical 
system  for  the  VARVS  is  to  quantify  the  angular  distortions  introduced  by 
the  1 ront  element.  From  tills  data  the  useable  aperture  can  be  determined 
which  in  turn  establishes  the  theoretical  resolution  limit  of  the  narrow 
field  systems.  To  accomplish  this,  a paraxial  meridional  ray  fan  covering 
the  entire  lens  clear  aperture  was  traced  through  the  front  element  of  the 
IR  VARVS  lens. 


The  angular  deviation  of  the  ray  fan  is  plotted  in  Figure  46  and  are 
denoted  by  circles.  Note  that  the  optical  power  is  convergent  out  to  about 
7 inches  of  clear  aperture.  Theoretically  this  should  be  correctable. 
However,  the  reversal  in  ray  angle  that  occurs  at  about  5 inches  clear 
aperture  diameter  would  probably  create  serious  optical  aberrations.  For 
this  reason  we  chose  a maximum  aperture  of  5.5  inch.  In  order  to  gain 
some  insight  into  the  practicability  of  correcting  the  deviations  shown 
in  Figure  46,  the  angular  deviation  required  at  the  corrector  surface  to 
produce  perfect  focus  at  the  desired  focal  length  were  computed  and  are 
plotted  in  Figure  46  against  the  aperture  diameter  where  the  rays  originate. 
These  are  denoted  by  the  triangles.  These  data  show  the  corrector  shape  to 
be  easily  realizable  over  the  desired  aperture.  This  aperture  theoretically 
established  a diffraction  limit  on  resolution  of  about  0.1  milliradian. 

The  desired  resolution  is  5 times  better  than  the  wide  field  system.  In 
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Figure  46  Optical  characteristics  of  the  front  element  for  paraxial  rays 
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Section  V the  wide  field  system  yielded  a threshold  resolution  of  about 
1 cycle/milliradian  which  represents  a discrete  angular  resolution  of 
approximately  0.6  milliradian.  The  desired  narrow  field  resolution  is 
therefore  about  0.12  milliradian.  This  value  is  compatible  with  the 
aperture  diameter  limit  discussed  above.  Using  0.12  milliradian  as  a 
resolution  goal,  the  f ield-of-vlew  and  focal  length  requirements  can  be 
calculated  if  image  size  is  defined.  In  order  to  use  a common  scanner 
for  both  wide  and  narrow  f ields-of-view,  both  Image  sizes  must  be  identical. 
This  size  has  been  established  in  the  study  of  Reference  2 is  0.72  inches. 
The  f ield-of-view  and  focal  length  are  calculated  as  follows:  If  0.12 
milliradian  resolution  is  desired,  the  detectors  must  have  an  angular 
separation(o)  of  less  than  about 


If  488  effective 


Resolution  Desired 
a Kell  Factor 

scan  lines  are  available, 


0.12 

1.4 


0.086  mr 


the  f ield-of-view  (6)  will  be 


(28) 


0 = .086  x 488  = 42  mr 
0 = 2.4° 

The  focal  length  (f)  is  therefore 


f - 


Image  Size  _ 0.72 
Field-of-view  0.042 


17.14  inches 


(29) 


In  order  to  allow  some  margin  for  error,  this  was  increased  to  20  inches 
which  yields  a 2°  field-of-view. 

Another  requirement  for  the  common  aperture  system  is  that  the  optical 
output  be  compatible  with  common  modular  scanners.  Again  utilizing  experience 
gained  in  the  wide  field  feasibility  study  (Reference  2),  space  must  be 
provided  identical  to  that  optical  system  and  is  as  shown  in  Figure  45. 

2.  OPTICAL  DESIGN 

Having  looked  at  the  wavefront  exiting  from  the  first  element,  the 
initial  design  approach  was  centered  on  the  concept  of  correcting  the 
errors  Introduced  by  the  first  element  with  another  spline  element  (called 
a corrector)  which  would  collimate  the  energy,  and  then  to  focus  this 
recollimated  wavefront  with  a classical  NFOV  IR  design.  Initially  a 20 
inch  focal  length  parabola  was  used  for  the  latter.  Some  work  on  this 
approach  made  it  evident  that  the  corrector  power  required  to  recollimate 
the  wavefront  was  a detriment  to  performance.  The  paraxial  focal  length 
of  the  first  element  is  21.68  inches.  Because  this  power  is  similar  to 
that  expected  of  the  final  design,  the  initial  approach  using  a parabolic 
reflector  was  discontinued,  the  reflector  was  removed  from  the  lens 
prescription  and  optimization  continued  using  the  first  element  for  power 
and  the  corrector  as  a zero  power  component  functioning  to  "clean-up"  the 
wavefront. 
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While  working  with  this  two  element  configuration  and  optimizing  the 
corrector,  relationships  began  to  be  established  between  angular  resolution, 
aperture  diameter  and  corrector  complexity  (i.e.,  number  of  splines  per 
surface).  These  relationships  are  shown  in  Figure  47.  These  data  show 
on-axis  resolution  of  0.16  milliradians  was  obtainable  out  to  a 3.5  inch 
aperture.  At  5 inches  of  aperture,  resolution  degraded  to  0.47  milliradians. 
Figure  46  shows  a rapid  change  to  front  element  angular  error  in  the  4-5 
inch  aperture  region.  Therefore,  it  is  very  possible  that  not  enough 
splines  were  used  in  this  region.  Five  splines  per  surface  are  the  maximum 
available  in  our  design  program. 
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Figure  47  Optical  performance  of  early  design  long  focal  length  lens 


We  verified  this  by  using  another  commercial  optical  design  program 
(ACCOS  V)  which  is  not  limited  to  number  of  splines  per  surface.  This 
program  was  able  to  correct  a 5.5  inch  aperture.  However,  the  program 
had  continuity  control  problems  when  operating  in  the  optimization  mode. 

It  only  worked  properly  for  rays  passing  through  the  center  portion  of 
the  splines  and  not  the  end  points.  For  this  reason  we  abandoned  ACCOS  V 
for  spline  corrector  element  design.  Figure  48  shows  a significant 
resolution  drop  with  field  angle.  This  indicates  that  a more  complex 
optical  design  will  be  required.  To  study  this,  lens  elements  were  added 
after  the  corrector  and  the  f ield-of-view  increased  to  2 degrees  for 
optimization  of  off-axis  imagery.  The  startine  design  after  reoptimizing 
the  axial  blur  is  noted  as  RN11  on  Figure  48.  The  angular  blur  or  resolu- 
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Cion  of  this  configuration  operating  with  a 5.5  inch  aperture  is  shown 
in  Figure  48  for  zero  and  1°  field  angles.  The  points  designated  as  RN19 
are  the  result  of  several  design  runs  for  optimization  of  the  off-axis  blur. 
RN26  is  the  result  of  further  work  with  the  prescription  directed  at 
reduction  of  the  size  of  the  central  obscuration. 

To  gain  more  correction,  lenses  were  added  in  the  region  of  the  image 
plane,  resulting  in  the  design  which  has  produced  the  greatest  overall 
resolution  (RN27).  However,  this  design  as  it  now  stands,  requires  an 
increase  in  track  length  between  corrector  and  image  to  provide  adequate 
room  for  folding.  Also  the  focal  length  of  this  design  is  10.57  inches 
which  is  too  short  for  a f ield-of-view  of  2°.  Further  work  is  being  done 
to  increase  the  track  length  and  the  focal  length.  The  aperture  ratio  will 
be  approximately  4.0  when  the  desired  focal  length  is  obtained. 

The  ACCOS  V lens  design  program  was  used  in  a parallel  optimization 
effort.  The  first  element  (common  to  the  wide  angle  system)  is  entered 
in  the  prescription  as  a spline  and  is  not  a variable.  As  long  as  the 
spline  is  not  varied,  ACCOS  V does  not  have  the  continuity  control  problems 
mentioned  earlier.  The  best  design  from  ACCOS  V as  of  this  report  is 
SPL  25.  This  design  is  the  narrow  f ield-of-view  system  shown  in  Figure  45. 
All  eight  surfaces  in  the  relay  system  are  aspheric,  four  carrying  deform- 
ation coefficients.  The  performance  data  on  this  design  for  comparison 
with  the  spline  corrector  system  is  plotted  in  Figure  48.  This  design  at 


Figure  48  Comparison  of  optical  designs 
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present  is  monochromatic  at  12  microns.  Spectral  correction  is  not  expected 
to  be  a problem  based  on  the  work  with  the  spline  program.  The  central 
obscuration  diameter  of  this  design  is  2.5  inches  and  is  defined  by  the 
lens  element  nearest  the  folding  mirror.  The  aperture  stop  needs  to  be 
moved  inside  the  first  element  and  located  at  the  first  corrector  surface. 

By  plotting  the  data  of  Figure  47  for  the  same  aperture  (RN8A)  on 
Figure  48  it  is  clear  that  optimizing  for  the  1°  field  angle  has  reduced 
on-axis  performance  well  below  our  requirements.  Furthermore,  any  optimiz- 
ation effort  aimed  at  improving  off-axis  performance  causes  a corresponding 
reduction  of  on-axis  performance  and  vice  versa.  For  this  reason,  the 
design  work  was  terminated  and  an  investigation  was  begun  to  find  the  reason 
for  the  low  performance.  The  source  of  the  problem  was  found  to  be  the 
wavefront  aberrations  introduced  by  the  spline  surfaces  of  the  first  element. 
The  corrector  can  accurately  correct  these  for  one  object  angle  but  some 
compromise  must  be  made  for  other  object  angles.  It  is  clear  that  either 
aperture  diameter  or  field  angle  must  be  reduced  to  obtain  satisfactory 
performance.  To  study  this  the  angular  relationships  introduced  by  the 
first  element  were  determined  for  on-axis  and  off-axis  conditions.  This 
was  accomplished  by  tracing  collimated  meridional  ray  fans  each  containing 
21  rays  through  the  element  at  0,  1/4°  and  1°  object  field  angles  as  shown 
for  greatly  exaggerated  angles  on  Figure  49.  These  data  were  used  to 
compute  the  amount  an  off-axis  ray  is  deviated  relative  to  the  same  ray 
on-axis.  For  any  ray  n: 

6 = A 6'  - a" 

n n n 

This  value  is  plotted  against  ray  position  in  terms  of  aperture  diameter  in 
Figure  50  for  the  two  object  angles.  We  were  very  surprised  to  find  that 
the  data  fell  along  a straight  line  for  both  field  angles  out  to  4.8  inch 
aperture. 


52 


Had  the  points  deviated  from  a smooth  curve  by  an  amount  exceeding 
the  system  angular  resolution  requirement  corrected  for  the  effects  of 
the  dovm-stream  optical  elements  required  to  give  the  correct  focal  length, 
optical  design  would  be  clearly  impossible.  The  adjusted  tolerance  is 
shown  on  Figure  50.  All  points  fall  within  the  tolerance  out  to  4.8  inch 
aperture.  While  the  linear  buildup  in  ray  deviation  with  aperture  is 
encouraging,  this  does  not  necessarily  mean  design  is  possible. 

The  slope  of  the  curves  on  Figure  50  are  an  indication  of  the  degree 
of  aberration  produced  by  the  optical  element.  An  aberration  free  2s 
f ield-of-view  optical  element  would  produce  the  nearly  vertical  dashed  line 
on  Figure  50.  The  difference  between  this  and  our  lens  is  very  large. 

This  indicates  very  large  aberrations  and  is  probably  the  source  of  our 
design  difficulty.  Nevertheless  we  feel  the  linear  build-up  in  aberrations 
with  aperture  diameter  should  make  them  correctable.  As  an  example,  rays 
through  a meniscus  lens  with  the  same  focal  length  as  our  lens  are  shown 
on  Figure  50  for  a 1°  field  angle.  This  curve  compensates  with  conventional 
spherical  optics. 

However,  because  of  time  constraints  we  had  to  terminate  the  effort 
at  this  point  and  Intend  to  continue  this  independently. 


Angular  Difference,  t>n  - deg 
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Figure  60  Optical  errors  for  0.26°  and  1 .0°  rays 
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3.  SUMMARY 


As  of  this  report,  we  have  not  been  able  to  obtain  an  optimized  design 
for  the  narrow  f ield-of-vie»  case  which  meets  the  resolution  goals.  How- 
ever, considerable  insight  has  been  gained  into  the  design  problems  and 
we  are  reasonably  confident  that  a narrow  field  optical  system  can  be 
added  to  the  VARVS.  However,  aperture  constraints  will  probably  limit 
its  resolution  to  approximately  0.2  milliradians. 


APPENDIX  A 


APPLICATION  OF  THE  NIGHT  VISION  LABORATORY  (NVL) 

THERMAL  VIEWING  SYSTEM  STATIC  PERFORMANCE  MODEL  TO 
THE  VISUAL  SPECTRUM  VARVS 

It  was  suggested  that  the  NVL  Thermal  Viewing  System  Static  Performance  Model 
be  used  to  evaluate  the  performance  of  the  Remote  Viewing  System  (RVS).  However, 
repeated  attempts  to  convert  the  RVS  parameters  directly  to  the  NVL  model  have  led 
to  the  following  problem.  The  radial  distortion  function  of  the  foveal  lens  does 
not  lend  itself  to  an  MTF  analysis  as  a function  of  object  field  angular  spatial 
frequency  as  called  for  in  the  NVL  model.  All  parameters  can  be  converted 
successfully  except  for  the  scan  velocity  term  because  a linear  raster  scan  on 
the  lens  image  plane  will  create  a variable  angular  velocity  and  variable  direction 
scan  in  the  object  field.  This  is  depicted  in  Figure  A-l.  Extreme  complexity 
results  when  attempts  are  made  to  convert  spatial  into  temporal  frequency.  This 
is  illustrated  by  the  rotation  of  the  f bar  pattern  in  the  lens  image  plane 
shown  in  Figure  A-l.  Given  enough  time,  an  analysis  could  be  made  in  a manner 
compatible  with  the  NVL  model.  However,  the  analysis  is  much  simpler  if  performed, 
not  in  object  field  angular  frequency  (cycles/milliradian)  but  in  spatial  frequency 
terms  (lines/millimeter).  For  our  purpose  of  optimizing  the  RVS  lens,  it  is 
simpler  to  work  in  terms  of  spatial  frequency  on  the  foveal  lens  focal  plane. 

This  simplicity  arises  because  seven  of  the  nine  MTF's  are  independent  of 
object  field  angle  at  this  foveal  lens  focal  plane  location,  and  the  scan 
velocity  is  undirectional  and  uniform  at  this  location,  thereby  making  easy 
conversions  from  spatial  to  temporal  parameters.  The  only  non-linear  conversions 
necessary  are  simple  geometrical  ones  which  translate  from  focal  plane  to  object 
field  and  display  space.  The  advantages  of  working  in  the  spatial  frequency 
terms  will  become  clear  as  the  analysis  is  developed.  In  the  following  develop- 
ment, the  NVL  model  approach  will  be  used  precisely  but  will  be  applied  in  the 
foveal  lens  focal  plane  as  a function  of  linear  spatial  frequency  (1/mm).  Parameters 
will  be  covered  in  the  same  order  as  they  are  in  the  NVL  Report  (4),  which, 
describes  the  model  in  detail. 


1.  MTF'S 
a.  Optical  MTF 

The  optical  MTF's  consist  of  a diffraction  MTF  and  a Gaussian  MTF. 

(1)  Diffraction 

In  angular  terms,  the  diffraction  MTF  is  reference  as  Equations  (9) 
and  (10)  of  the  NVL  report: 

H _(f  ,6)  - - [cos"1  A - A(1  - A2)1/21  (A-l) 

opt'  x*  it  L J 

where  A • XF^fx/L(0)  (A-2) 
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Figure  A-1  Scan  distortion  introduced  for  foveal  lens 


GP7t  103M0* 


where  L(0)  is  the  equivalent  focal  length  which  changes  over  a 50/1  range  as 
object  field  angle  6 changes.  The  angle  0 is  the  absolute  angle  between  the 
point  of  interest  and  the  lens  optical  axis.  At  the  foveal  lens  image  plane 

f 


where  S is  the  image  plane  spatial  frequency  and  f is  its  object  field 
angularXequivalent  measured  along  the  scan  line  projection  in  the  object 
field  (p  direction  on  Figure  A-l).  Solving  for  f in  Equation (A-3j  and 
substituting  this  for  f^  in  Equation  (A-2).  U 

A ■ IF,  S 


(A-3) 


Since  the  F/number  of  our  lens  is  constant,  the  diffraction  MTF  is  no  longer  a 
function  of  object  field  angle.  Thus  we  may  write  H t(S  ) which  indicates 
that  the  MTF  is  a function  of  the  independent  variab?lCS  xonly.  Note,  however, 
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that  conversion  to  object  field  angular  spatial  frequency  is  very  simple 
because  focal  length  is  constant  over  small  angular  increments  and  may  be 
determined  from 


fw  “ SxL(6)  (A- 5) 

where  p is  along  the  scan  line  projection  in  the  object  field 
likewise 

fw  = SyL(0)  (A-6) 

where  w is  normal  to  the  scan  direction  in  the  object  field 

(2)  Blur 

A similar  simplicity  exists  here.  The  MTF  equation  with  the  angular 
term  b of  Equation  (11)  of  Reference  ( 4 ) replaced  with  its  equivalent  is: 

wv>-~p[-Sf‘2]  (A‘7> 

4 

The  foveal  lens  inherently  has  a constant  spatial  blur  over  its  entire  focal  plane, 
so  that  the  sigma  (a)  of  Equation (A-7) is  a constant.  Substituting  Equation(A-5) into 
(A-7)  we  see  the  blur.  MTF  simplifies  to 

Hblur(Sx)  “ exP  C-  2"2  <A“8) 

Thus  this  MTF  like  the  diffraction  MTF,  is  no  longer  a function  of  object  field 
angle  because  the  focal  length  variable  has  been  removed. 

b.  Detection  MTF 

The  spatial  filter  MTF  of  the  detector  is  defined  as: 


Sin (nf  Ax) 

■ “» et ■ — ;-5-  ■ ““‘V*  <A-9) 

It  is  also  complex  in  our  system  because  the  angular  projection  of  the  detector 
into  the  object  field  (Ax)  in  this  equation  varies  with  absolute  object  field 
angle  (0).  Since  the  detector  height  is  still  uniform  at  the  lens  focal  plane, 
shown  in  Figure(A-2) as  Ah,  Equation(A-9)can  be  restated  as: 


Sin(itS  Ah  ) 
HDet(Sx)  “ vrA-h;— 


(A— 10) 
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Again  the  MTF  becomes  independent  of  object  field  angle.  Note  from  Figure (A-2)that 
the  detector  height  (Ah^)  is  a function  of  detector  size(a),  detector  system  focal 

length  (Lp),  and  relay  focal  length  (Lc),  viz: 


Ahx  * (A-ll) 


If  the  detector  characteristics  are  known,  the  focal  lengths  are  a function 
of  detector  size  (Ah)  projected  unto  the  image  plane  as  shown  in  Figure  (A-2). 
Detector  size  Ah  can  be  computed  directly  from  either  the  on-axis  resolution 
required,  the  number  of  scan  lines  required  across  the  vertical  FOV,  or  bandwidth/ 
response  restrictions  and  frame  rate  requirements.  The  focal  lengths,  L and  Ln, 

are  then  selected  to  make  the  detector  dimension  appear  as  the  required  Ah  at 
the  foveal  lens  focal  plane. 


MTF  becomes : 


WV  • S1"‘  rr 

x D 


Again  this  MIF  is  Independent  of  object  field  angle. 


(A-12) 
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c.  Detector  Electronics  MTF 


It  is  in  the  MTF,  the  detector  electrical  response,  that  we  get  into 
real  trouble  trying  to  work  in  object  field  angular  space.  For  a conven- 
tional linear  optical  system,  a linear  detector  scan  velofcity  converts 
into  a scaled  but  linear  angular  scan  in  the  object  field.  This  is  not 
true  in  our  system  as  was  shown  in  Figure  A-l.  A linear  scan  in  the  x 
direction  on  the  image  plane  results  in  angular  velocities  in  both  6X  and 
6V  directions  in  the  angular  object  field.  Both  of  these  angular  components 
are  nonlinear  functions  of  both  x and  y position  on  the  image  plane.  Thus, 
converting  from  spatial  frequency  to  temporal  frequency  becomes  very 
complex.  All  of  this  can  be  avoided  by  working  in  linear  spatial  plane 
terms.  If  the  scanner  has  an  angular  scan  velocity  B,  then  the  linear 
motion  of  the  instantaneous  FOV  on  the  foveal  lens  image  is 

Vx  - BLC  (A-13) 

The  conversion  to  temporal  frequency  (f)  is  therefore 

f = V S (A-14) 

xx 

This  is  a constant  conversion  and  not  a function  of  time.  Therefore,  all 
electronic  MTF's  of  the  NVL  model  are  valid.  These  are 


H’Det<f) 

HElect(f) 


Vf> 


d.  Display 

The  RVS  display  is  the  inverse  of  the  foveal  lens,  which  results  in 
a conventional  linear  raster  generated  on  the  CRT.  The  CRT  has  a constant 
spot  size  and  the  expansion  optics  has  a constant  blur  at  the  object  focal 
plane.  Again  this  MTF,  if  derived  in  the  linear  spatial  plane,  will  not  be 
a function  of  object  angle.  If  the  optical  blur  and  CRT  spot  size  are 
combined  and  assumed  to  have  a Gaussian  MTF,  a composite  sigma  (oj)  results 
and  the  MTF  is: 


HDisp(Sx)  “ eXp  [-2*2<r°d)2sx2  J (A"15) 

where  r is  the  physical  ratio  of  format  sizes  ; viz 

^LENS  IMAGE  (A-16) 

|-  m i ■ . ■ x 

Display  crt 
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By  contrast,  if  this  were  accomplished  in  the  object  angular  plane,  the  MTF 
would  be  much  more  complex,  viz 


HDisp(fx,0,M) 


2"2<”d)2fx21 

2 2 J 

ne)V  J 


(A-17) 


where  M is  any  system  angular  magnification  from  object  field  to  the  viewer. 
Again  the  simplicity  is  obvious. 

e.  & f.  Stabilization  and  Eyeball 


The  remaining  two  MTF’s  are  the  only  two  that  are  not  aimplified  by 
working  in  linear  spatial  rather  than  angular  terms.  First,  stabilization 
tends  to  be  angular  input  to  the  system.  Using  the  MTF  from  the  NVL  report: 


V.(tx>  ■ “P  '-pfx2) 

Converting  to  the  foveal  lens  image  plane  results  in 

HLos(Sx,G)  ‘ exp  C-PS*2l<G>2]  ‘ (A-19) 

Similarly,  the  eye  views  the  display  in  angular  terms;'  The  NVL  MTF  is 


HEye(fx)  " exp["  1 r] 

Equation (A-20)must  be  converted  to  the  foveal  lens  image  plane 


(A-20) 


(A— 21) 


In  conclusion,  seven  MTF's  have  been  simplified  at  the  expense  of  two  that 
have  been  made  slightly  more  complex  by  the  conversion  to  linear  spatial 
frequency. 

2.  NOISE  EQUIVALENT  MODULATION  (NEM) 

For  visual  spectrum  applications  noise  equivalent  modulation  must  replace 
NE A T in  the  NVL  model.  In  the  visual  model,  the  primary  noise  source  is  the 
detector  which  is  a silicon  vidicon.  Its  NEM  was  extracted  from  data  of 
Reference  ( 7),  These  data  show  vidicon  S/N  as  a function  of  faceplate  illumin- 
ation for  a specific  bandwidth.  The  basic  function  is  approximately 

Peak-to-peak  signal  - 100  E (A-22) 

noise  (rms) 

where  E is  faceplate  illumination  in  LUX.  The  noise  equivalent  signal  is 
(signal  input  that  just  equal  noise) 


NEM  - 


noise 

signal 


1 

100E 


(A-23) 


assuming  that  the  noise  is  proportional  to  the  square  root  of  the  bandwidth 
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f)  of  4(10®)  Hz.  For  data  given: 

A f 


NEM 


6 - 5 X 10-6  (E  in  LUX)  (A-24) 


For  E in  footcandles: 


NEM  - 


100E  4X  10 


4.64  x 10" 7 A f (E  in  Foot-Candles)  (A-25) 


The  faceplate  illumination  can  be  calculated  from  system  geometry  as  follows: 

(A-26) 


Ef  - 


_T  T 

B a o 


4 F4 


No 


Where 


B - Scene  brightness  in  footlamberts 

T ■ Atmospheric  transmission 

■ Optical  transmission  with  sensor 

F ■ The  equivalent  F/number  or  F/number  actually  supplying  the 

vidicon.  This  is  the  lens  F/number  modified  by  the  relay  and 
from  basic  geometrical  optical  theory  is: 


F 

noe 


Fno  LD 


(A— 27) 


If  the  sensor  employs  an  automatic  light  level  control  which  operates  on 
vidicon  target  current,  E will  be  accurately  maintained.  Therefore,  Equation 
(A-25)  applies  as  written  for  the  level  of  E which  is  preset.  For  the  silicon 
vidon  under  study,  best  performance  is  obtained  when  the  level  is  about 
2 

0.1  lumens/ft  . Equation  (A-24)  then  becomes: 

NEM  - 4.64  X 10"6  Af  (A-28) 


f 


I 


J 
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3.  MRT  CALCULATIONS 

Tl.e  following  MRT  equation  modifications  are  required  so  that  the  computation 

may  be  performed  in  linear  spatial  frequency  terms.  First,  in  the  NVL  MRT  equation,  Ay 

must  be  replaced  by  the  apparent  detector  sire  at  the  foveal  lens  image  plane, 

i.e.,  it  must  be  the  Ah  defined  on  Figure <A-2).  As  previously  demonstrated  in 
equation  (A—ll). 

L„ 


Ah 


(A-29) 


Also,  in  the  MRT  equation,  it  is  best  to  compute  the  Q integral  in  terms  of 
temporal  frequency.  This  eliminates  the  velocity  term  in  the  MRT  equation 
and  makes  the  Q integral  easier  to  compute.  The  Q integral  is  therefore 


Q(f 


*8)  ~ J 


(A-30) 


Of  these  terms,  only  H^,  the  transfer  function  for  a rectangular  bar  of  width  w, 

has  not  been  defined.  This  transfer  function  is  in  linear  rather  than  angular 
dimensions,  i.e.. 


«„(!*)-  Sl»c(w£x)  - Sl„c(usx) 

' V / v 


(A-31) 


where 


u A 1 

W m — — 

2S 


(A-32) 


Ihe  MRT  equation  written  to  show  the  dependency  of  two  variables  is 


MRT(Sx,8) 


SNRtrNEAT 

4^14  MTF_„„  , (s  ,0) 
TOTAL ' x ' 


Ahy  SxQ(f,6) 
,^fNFR  te  nOVSC 


1/2 

(A-33) 


Ihis  equation  results  in  an  MRT  very  weakly  dependent  on  6.  To  obtain  the  MRT 
for  any  field  angle  e we  convert  the  spatial  frequency  term  Sx  into  an  angular 

frequency  term  by  using  Equation (A-5)containing  the  focal  length  function: 


{u  - SxL(6)  (A-5) 

Note  this  will  be  the  angular  spatial  frequency  in  the  scan  direction  (target 
bars  normal  to  the  scan  direction).  It  could  be  related  to  f and  f but 

X y 

this  does  not  appear  to  be  required  at  this  point. 
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4.  CONCLUSIONS 


To  conclude  this  effort,  a block  diagram  of  the  NVL  model  converted  to  the 
VARVS  Concept  in  the  visual  spectrum  is  shown  in  Figure  A-3.  This  model  was 
used  in  the  study  to  compute  Minimum  Resolvable  Modulation  to  predict  performance. 


Figure  A-3  NVL  model  adapted  to  VARVS  for  visual  spectrum 


APPENDIX  B 


LINEAR  SENSOR  ACQUISITION  MODEL 


I 

l 
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The  linear  sensor  acquisition  model  developmer t was  performed  with 
IRAD  funds  and  is  included  in  this  report  for  reference.  The  reader 
will  find  the  model  framework  similar  to  that  of  the  VARVS  search  model 
The  main  differences  lie  in  the  linear  sensor's  relatively  constant 
resolution  over  the  f ield-of-view  thus  enabling  the  observer  to  make 
better  use  of  his  time  by  eliminating  the  need  to  slew  the  sensor,  for 
target  identification. 


This  model  was  designed  to  evaluate  the  use  of  TV  sensors  for  target 
acquisition  and  computes  the  probability  of  cue  detection,  target  detection, 
recognition,  and  probability  that  the  target  is  within  the  search  area.  The 
model  requires  input  parameters  defining  the  mission  geometry,  search  area 
requirements  and  sensor  and  sensor  system  characteristics.  The  model 
assumes  a cue  oriented  systematic  search  which  is  one  where  targets  are 
found  by  first  detecting  a cue  and  then  using  the  cue  to  search  for  a 
target.  For  example,  a tank  or  personnel  carrier  may  be  located  by 
finding  a cluster  of  vehicles.  The  probability  of  detecting  a cue  or  target 
has  been  modeled  as  an  exponential  as  a function  of  resolution  elements 
across  the  critical  (smallest)  dimension  of  the  cue  or  target.  The  model 
allocates  the  time  available  for  search  so  that  the  combination  of  tasks 
performed  to  acquire  a target  is  optimized. 

The  main  parameters  controlling  the  search  process  are  derived  from  the 
TV  sensor's  footprint  which  is  a description  of  the  ground  resolution  of 
the  sensor  referenced  to  the  detector  plane.  The  resolution  measurement 
has  units  of  line  pairs  per  millimeter.  Characterization  of  the  footprint 
is  accomplished  by  calculating  system  resolution  at  uniform  increments 
(AA)  on  the  detector  face.  An  adaptation  of  the  NVL  model  is  used  for 
this  calculation.  The  quality  of  the  footprint  affects  the  methods  and  use 
of  the  footprint  during  search.  The  resolution  will  vary  throughout  the 
footprint  which  is  a result  of:  1)  the  slant  range  to  each  AA,  2)  the 
image  degradation  due  to  the  forward  velocity  of  the  aircraft,  3)  atmosphere, 
and  4)  viewing  geometry. 


The  user  may  determine  that  parts  of  the  display  are  unsuitable  for 
his  search.  Examples  of  this  include  information  displayed  above  the 
horizon,  ground  covered  by  sensor  which  has  too  long  of  a slant  range  so 
that  target  resolution  is  below  the  threshold,  and  other  imagery  which 
is  below  the  resolution  threshold  due  to  image  motion  or  masking  problems. 
Thus,  the  user  may  decide  to  ignore  portions  of  the  display  for  reasons 
stated  above  or  other  reasons  which  become  significant  in  particular 
applications.  The  net  effect  of  ignoring  parts  of  the  display  is  to 
limit  the  "real"  coverage  of  the  sensor. 

The  size  and  shape  of  the  real  coverage  of  the  sensor  plays  a major 
part  in  its  dynamic  search  capabilities.  As  the  size  of  the  sensor's 
useable  footprint  is  decreased,  the  demands  of  the  dynamic  characteristics 
of  the  sensor '8  servo  and  search  process  is  Increased.  This  suggests  that 
a tradeoff  exists  for  the  optimum  choice  of  the  portion  of  the  footprint 
used.  This  choice  may  vary  with  the  specific  mission  parameters.  An 
examination  of  the  dynamic  systematic  search  process  will  give  some  insight 
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into  the  necessary  parameters  to  consider  in  this  tradeoff.  A typical 
linear  sensor  is  shown  in  Figure  B-l.  The  outer  edges  of  the  footprint 
are  determined  by  the  viewing  geometry  and  sensor  FOV.  If  the  user 
determines  that  the  lower  portion  of  the  footprint  is  to  be  used  in  search, 
the  real  coverage  can  be  characterized  by  the  length  along  the  flight  line 
Xp  and  the  mean  width  2Yp.  In  a systematic  search,  the  sensor  will  typically 
be  slewed  in  a snapshot  fashion  across  the  total  search  area.  The  term 
"snapshot"  indicates  that  the  sensor  is  slewed  quickly  to  its  destination 
and  then  held  stationary  while  the  user  searches  the  display.  This  type 
of  slewing  is  used  to  eliminate  smearing  on  the  detector  due  to  the  angular 
velocity  of  the  sensor. 


S«arch 

Geometry 


Typical  Ground 
TV  Footprint 


20°  FOV  TV 
2000  ft  Altitude 
6000  ft  Downrange 


3500  ft 


Figure  B-1  Linear  sensor  search  geometry 


To  assure  full  coverage,  the  footprints  must  fill  the  required  lateral 
coverage  2YT  during  the  time  the  sensor  platform  moves  X_  along  the  flight 
line  is  illustrated  in  Figure  B-2.  The  number  of  slews  (n8i)  required  to 
cover  the  required  offset  2Yf  is  defined  as: 


n 


si 


(B-l) 
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where  []  indicates  rounding  up  the  next  whole  number.  The  time  available 
(ta)  to  find  a particular  target  is  calculated  in  terms  of  the  sensor  plat- 
form velocity  (V^) , 


X 


The  time  it  takes  the  sensor  to  move  to  the  new  position  and  settle  to  its 
steady  state  value  will  be  called  sensor  slew  time  (tsi) . The  amount  of  time 
the  observer  will  spend  "searching"  during  a particular  "snapshot"  is  related 
to  his  glimpse  rate.  The  time  per  glimpse  (tg)  has  been  measured  to  be 
between  0.25  to  0.3  seconds. 


The  time  (ta)  required  for  a complete  swath  of  the  acquisition  process 
is: 


t' 

a 


si  v si 


+ n 


V 


(B-3) 


where  ng  is  the  number  of  glimpses  per  slew  and  n8^  is  the  number  of  slews. 

If  the  user  is  required  to  keep  a constant  nominal  depression  angle  it 
is  apparent  that 


If  t'  > t , the  user  must  be  "tracking"  the  search  area.  This  restriction 
brings  to  light  a second  tradeoff  or  optimization  that  must  be  made.  The 
user  may  elect  to  spend  more  time  searching  each  "snapshot"  and  not  cover 
the  total  search  width.  This  strategy  becomes  important  when  the  total 
slew  time  (n8i  • ts^)  becomes  significant  with  respect  to  the  pure  search 
time,  (ngl  • ng  • tg). 
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The  time  allocation  tradeoff  and  the  footprint  coverage  compromise 
constitutes  the  basic  mechanics  of  the  model.  The  performance  measure  which 
determines  how  the  two  tradeoffs  will  be  made  is  a function  of  conditional 
probability  functions.  These  functions  define  the  probabilistic  description 
of  the  particular  tasks  in  terms  of  the  number  of  resolution  elements 
across  the  critical  dimension  of  the  cue  or  target. 

The  probability  functions  are  exponential  of  the  general  form 

P(y)  - 1 - exp  [-(a  y + b)2]  (B-4) 

where  y = no.  of  resolution  elements  across  the  target 

a,b  = constants 

Three  tasks  are  defined  by  these  functions.  They  are: 

1.  The  probability  of  detecting  the  cue  given  that  you  look  at  the 
cue,  (P'), 

2.  The  probability  of  detecting  the  target  given  that  you  have 
detected  the  cue,  (P^,) , 

3.  The  probability  of  recognizing  the  target  given  you  have  detected 
the  target  (Pr). 

The  first  task  is  a function  of  Nc,  the  number  of  resolution  elements 
across  the  critical  dimension  of  the  cue.  It  is  generally  accepted  that  one 
line  pair  across  the  target  produce  a P'(l)  equal  to  0.90.  Also  no 
information  is  available  to  the  observer  when  the  available  resolution  is 
1/7  to  1/8  elements  across  the  cue,  Pc(0.13)  * 0.  These  two  values  can  be 
used  to  define  a and  b in  Equation  B-4.  Then  the  conditional  probability 
function  associated  with  Task  One  is  defined  as: 


P'(N  ) = 1 - exp  [-(1.75  • N - ,23)2]  (B-5) 

c c c 

Task  two  has  a similar  derivation  and  the  only  difference  is  that  Nt, 
the  number  of  resolution  elements  across  the  target  is  used  for  its 
argument.  Thus: 

PT(Nt)  - 1 - exp  [-(1.75  Nt  - .23)2]  (B-6) 

The  probabilistic  description  of  Task  Three  is  taken  directly  from  Bailey's 
paper  ( 6 ) and  is  defined: 

Pr(Nfc)  - 1 - exp  [— (*2 lJ  1 (B-7) 

Equation  B-5  defines  the  conditional  probability  of  cue  detection 
for  one  cue  and  one  glimpse,  and  does  not  take  into  account  the  possibility 
of  multiple  glimpses. 
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The  effect  of  multiple  glimpses  and  multiple  targets  was  studied  in 
detail  before  deciding  on  the  type  of  model  to  be  used.  The  performance 
of  the  observer  can  best  be  analyzed  by  answering  the  question,  "What  is  the 
probability  of  finding  exactly  K targets  in  ng  glimpses?".  This  question 
is  usually  associated  with  Bernoulli  Trials. 


P(K)  = (J)  pK  (1  - p)N_K 

where  <B-8> 

A = N! 

V (N-K)!  K! 

In  the  multiple  target  acquisition  case  the  probability  of  finding  a cue 
is  a function  of  the  number  of  unfound  cues  in  the  display.  Therefore  as 
the  observer  detects  cues,  the  probability  of  finding  another  cue  changes. 

If  the  probability  of  finding  a cue  per  glimpse  is  (P£)  for  one  cue  contained 
within  the  display,  then  the  probability  (P(NT))  of  finding  a cue  where 
multiple  (NT)  cues  exist  is: 

P(NT)  = 1 - (1  - P')NT  (B-9) 

Where  NT  is  the  number  of  independent  cues  contained  within  the  display.  As 
the  observer  detects  and  checks  out  possible  cues,  P(NT)  decreases  because 
the  number  of  unfound  cues  is  decreasing.  Due  to  this  non-stationary 
characteristic  of  P(NT),  the  expected  number  of  cues  E(NT)  found  during  K 
glimpses  cannot  be  calculated  using  compact  methods  such  as  Bernoulli  trials 
(a  form  of  the  Binomial  expansion).  In  the  *target  acquisition  application 
considered  here,  the  probability  density  function  defines  the  probability  of 
finding  K cues  in  ng  glimpses  for  constant  ng  and  all  K such  that  K is  whole 
number  0 _<  K <_  ng.  Then  one  must  determine  all  possible  ways  to  detect  K 
cues  in  ng  glimpses  (all  permutations),  find  the  probability  associated  with 
each  permutation  and  sum  these  probabilities.  With  the  use  of  a digital 
computer  and  its  memory  capabilities  this  process  may  be  somewhat  simplified. 
Consider  the  state  diagram  in  Figure  B-3.  The  circles  represents  specific 
states  as  defined  by  the  two  axes,  number  of  glimpses  made  (trials)  and  cues 
detected  (successful  results). 

Each  state' has  associated  with  it  a probability  of  being  in  that  state 
given  the  number  of  glimpses  made.  The  probability  of  traveling  from  one 
state  to  another  state  is  denoted  by  the  flag  on  the  arrow  joining  the 
two  states.  Each  node  defines  the  discrete  probability  density  function 
of  finding  exactly  K cues  in  n^  glimpses. 

To  determine  the  probability  associated  with  any  specific  state, 

Pc(ng,K),  one  must  determine  all  possible  paths  from  (0,0)  to  the  state 
in  question.  The  probability  of  traveling  down  a specific  path  is  the 
product  of  all  the  flags  along  that  path.  The  summation  of  all  possible 
probabilities  associated  with  various  path  ways  defines  Pc(ng,K).  The 
probability  of  being  in  state  (0,0)  is  P(0,0)-1. 


, 
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GP7I-1  til-6 

Figure  B-3  Acquisition  probability  density  function 


The  calculations  required  to  generate  the  density  distribution  can  be 
simplified  if  one  makes  the  observation  that  each  state  Pc(ng,K)  is  a function 
of  at  most  two  previous  states,  for  example: 


Pc(3,2)  = Pc(2,2)  • p(NK  - 2)  + Pc(2,l)  - p(K-l)  (B-10) 


Pc(l,0)  - Pc(0,0)  • p(K) 

Using  a digital  computer  and  saving  the  appropriate  information,  the  trials 
matrix  can  be  built  up  with  a reasonable  number  of  calculations. 

The  expected  number  (E[K(ng)]  of  possible  cues  found  is  by  definition 


E[K(ng)l 


oo 

*/k 


P(ng)dK 


(B-ll) 


or  because  of  the  discrete  nature  of  this  function. 


NG 

E[K(ng)  ] " 1 ' P(ng,D 

i-1 


(B-12) 
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If  all  possible  cues  have  equal  probability  of  being  the  true  cue,  the 
probability  of  finding  the  right  cue  Pc  is  simply 

p = E[K(ng) ] (B-13) 

C NT 

The  number  of  glimpses  the  user  is  able  to  make  is  both  a function  of 
the  average  glimpse  rate,  the  sensor  slew  time  and  the  number  of  targets 
within  the  display  on  average,  as  shown  in  Equation  B-3.  The  maximum  number 
of  glimpses  the  observer  can  make  is: 

",  ■ (B-14) 

The  final  probability  (P8)  to  be  calculated  is  the  probability  that 
the  cue  (target)  was  at  the  same  point  in  the  search  within  the  field  of 
view.  Assuming  the  probability  of  target  location  is  uniform  over  the 
search  area,  P can  be  calculated: 

w c 


Area  covered 

Total  Search  Area 


(B-15) 


The  area  covered  (Ac)  in  an  efficient  systematic  search  is  proportional 
to  the  total  number  of  slews  (N^  ) made: 


A 


c 


(B-16) 


The  term  Afp  represents  the  area  of  the  footprint.  The  total  search 
area  (Ax)  is  assumed  to  be  a rectangle  of  length  Xx  along  flight  path  and 
width  2Yx*  The  total  area  is: 


AT  - 2 . YT  • XT  (B-17) 

The  number  of  slews  the  operator  is  able  to  make  during  the  total  search 
process  is  defined  in  terms  of  the  platform  velocity,  search  area  length, 
and  time  taken  per  "snapshot"  and  is: 


V 

x 


+ n 


V 


Substituting  into  Equation  (B-15): 


(B-18) 


PS  " 


JJL 


2 Vx  <‘sl  + ”g  V YT 


(B-19) 


Equation  (B-19)  shows  that  under  the  assumption  made,  systematic  search 
with  a constant  depression  angle,  the  length  of  the  search  area  has  no 
effect  on  the  optimum  search  performance.  In  a specific  mission  with 
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constant  ng  glimpses  per  slew,  as  Y^.  increases  P8  will  decrease.  This 
decrease  in  P8  may  be  offset  by  decreasing  ng  thus  increasing  P8,  and 
decreasing  the  calculated  value  of  Pc.  The  model  adjusts  these  terms  for 
best  results.  In  light  load  situations,  for  example  low  aircraft  velocity 
and  "fast"  slew  time,  the  upper  bound  of  Y-j  may  be  the  geometric  crossrange 
capability  of  the  sensor  due  to  depression  angle  specifications,  and  not 
restricted  by  operator  bandwidth. 


A overview  of  the  acquisition  model's  implementation  is  shown  in  the 
flow  chart  in  Figure  B-4. 


The  model's  first  task  is  to  generate  the  sensor  footprint.  The  user 
is  required  to  designate  the  sensor's  geometric  location,  velocity,  atmos- 
pheric conditions,  masking  parameters  and  sensor  f ield-of-view.  Parameters 
defining  the  internal  characteristics  of  the  sensor  are  internal  to  the 
model  and  are  adjusted  by  modifications  to  the  actual  program. 

The  next  step  is  to  define  the  portion  of  the  footprint  that  is  to  be 
used  in  search.  The  model  requests  from  the  user  the  number  of  resolution 
elements  across  the  cues  to  be  used  as  a threshold  for  the  augmented 
footprint  choice.  The  augmented  footprint  is  then  characterized  in  terms 
of  length,  width,  and  average  resolution.  These  characteristics  are  used 
to  calculate  the  time  available  for  search,  the  conditional  probability 
of  detecting  the  cue  in  a single  glimpse  (P').  The  allocation  of  time  for 
search  and  sensor  slewing  is  adjusted  for  optimum  search  performance. 

Upon  completion  of  this  task  the  model  calculates  the  conditional  prob- 
abilities of  target  detection  (Pt)  and  target  recognition  (Pr).  The 
probabilities  of  mission  success  are  evaluated  in  terms  of  the  probability 
of  detecting  the  target 

P » p • P,  . p (B-20) 

aet  c d s 

and  the  probability  of  recognizing  the  target 


P 


t 


P 

r 


( B— 2 1 ) 


The  probabilities  of  mission  success  are  outputs  for  the  users  evaluation. 

The  user  may  now  elect  to  use  these  results  of  search,  or  select  an 
alternate  definition  of  the  augmented  footprint.  The  program  calculates 
the  optimum  probabilities  of  target  detection  and  recognition  for  a given 
mission  geometry.  It  was  recognized  that  a great  deal  of  information  can 
be  reused  to  describe  mission  geometries  varying  only  in  depression  angle. 

This  concept  has  been  Implemented  in  the  linear  sensor  model.  The  model  is 
able  to  output  predicted  mission  results  for  depression  angles  in  degree 
increments  over  the  total  useful  range  of  the  sensor  for  very  little  increased 
computation  cost. 


To  do  this  the  actual  computer  model  is  broken  up  into  two  programs. 
The  first  generates  the  resolution  capabilities  for  the  indicated  sensors 
over  the  total  useful  sensor  range.  The  second  exercises  the  search  model 
and  optimizations  described. 
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Figure  B-4  Computer  flow  diegrem  for  lineer  tensor  terget  acquisition 
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APPENDIX  C 


THE  EFFECT  OF  NUMBER  OF  MEASUREMENTS  ON  SAMPLE 
PROBABILITY  ACCURACY 


The  accuracy  with  which  the  performance  of  an  experiment  can  be 
P’-i dieted  is  a function  of  the  number  of  samples  that  are  taken.  The 
greater  the  number  of  samples,  the  greater  is  the  confidence  of  the  computed 
probability.  An  estimate  of  the  number  of  samples  that  must  be  taken  can 
be  calculated  using  the  Chebyshev  Inequality.  The  basic  form  of  the 
inequality  is: 


o 

P [ | x | > e]  < [x  ] , for  any  e > 0 (C-l) 

2 

E 

In  other  words,  the  probability  (Pr)  that  the  absolute  value  of  the  random 
variable  x is  greater  than  or  equal  to  e is  less  than  or  equal  to  the 
expected  value  of  x^  divided  by  e2. 

In  the  experiment  the  outcome  of  a trial  (y^)  is  defined  as: 

y^  = / 1 if  the  target  is  found  (C-2) 

( 0 otherwise 

The  probability  of  these  two  actions  is  defined  in  terms  of  the  binominal 
distribution: 


and 


P(y±  = 1)  = P 

P(yi  =0)=q=l-p 


(C— 3 ) 


All  trials  (y^)  are  independent,  therefore  the  random  variable  x can  be 
defined  as  the  difference  between  the  random  variable  y^  and  its  mean 


(C-4) 


Substituting  into  Equation  (C-l)  results  in 


n 

n 

n 

% > i i £ 

[i  £’i- 

i-i 

i-1 

i-1 

i.  right  hand  side 

of  Equation  (C-5)  can  be 

expanded  as 

n m 

i 2 2 IyiyJ ' yi<y:)  ■ yj<yi)  + yi2’ 

i-1  J-l 


(C-6) 
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The  average  of  any  product  where  the  subscripts  are  not  the  same  is  zero 
because  all  the  y^’s  are  independent.  This  Equation  (C-6)  becomes 


j,  t«7 -.■>  !i: 

i-1  i-1 

Since  all  y^  are  identically  distributed: 


(C-7) 


u 

^ o2  - n • 
yi 


i-1 

The  variance  of  a binomial  distrubtion  is: 

2 

0 = p * q 

y 

Substituting  Equations  (C-7)  and  (C-8)  into  (C-5)  gives: 


Pr  [ |x|  > e]  < £-^-| 

n e 

If  we  define  the  left  side  of  Equation  (C-10)  as  y we  find  that 


(C-8) 


(C-9) 


(C-10) 


< p..-  q 


(C-ll) 


n e 


Then  by  rearranging  terms,  the  number  of  trials  n is  defined  as 


n > 


p :..q 

2 

y • e 


Pi1  - P) 

2 

y • e 


(C-12) 


Thus,  the  number  of  trials  (n)  is  defined  in  terms  of  the  probability  of  a 
successful  mission  (p),  the  tolerance  of  the  measurement  (e)  and  the 
probability  of  being  within  the  tolerance.  Equation  (C-ll)  is  graphed 
on  Figure  14  for  two  values  of  e,  the  measurement  tolerance.  An  inter- 
esting observation  is  that  the  number  of  trials  required  to  predict  Y increases 
to  a maximum  as  the  probability  of  Y (p)  tends  toward  0.5. 


75 


APPENDIX  D 


APPLICATION  OF  THE  NIGHT  VISION  LABORATORY  (NVL) 

THERMAL  VIEWING  SYSTEM  STATIC  PERFORMANCE  MODEL  TO 
THE  IR  VARVS 

It  was  suggested  that  the  NVL  Thermal  Viewing  System  Static  Performance 
Model  D-l  be  used  to  evaluate  the  performance  of  the  RVS.  However,  repeated 
attempts  to  convert  the  RVS  parameters  directly  to  the  NVL  model  have  led  to 
the  following  problem.  The  radial  distortion  function  of  the  foveal  lens 
does  not  lend  itself  to  an  MTF  analysis  as  a function  of  object  field  angular 
spatial  frequency  as  called  for  in  the  NVL  model.  All  parameters  can  be 
converted  successfully  except  for  the  scan  velocity  term  because  a linear 
raster  scan  on  the  lens  image  plane  will  create  a variable  angular  velocity 
and  variable  direction  scan  in  the  object  field.  This  is  depicted  in  Figure 
D-l.  Extreme  complexity  results  when  attempts  are  made  to  convert  spatial 
into  temporal  frequency.  This  is  illustrated  by  the  rotation  of  the  fx  bar 
pattern  in  the  lens  image  plane  shown  in  Figure  D-l.  Given  enough  time,  an 
analysis  could  be  made  in  a manner  compatible  with  the  NVL  model.  However, 
the  analysis  is  much  simpler  if  performed,  not  in  object  field  angular  fre- 
quency (cycles/milliradian)  but  in  spatial  frequency  terms  (lines/millimeter). 
For  our  purpose  of  optimizing  the  RVS  lens,  it  is  simpler  to  work  in  terms 
of  spatial  frequency  on  the  foveal  lens  focal  plane. 

This  simplicity  arises  because  seven  of  the  nine  MTF’s  are  independent 
of  object  field  angle  at  this  foveal  lens  focal  plane  location,  and  the 
scan  velocity  is  undirectional  and  uniform  at  this  location,  thereby  making 
easy  conversions  from  spatial  to  temporal  parameters.  The  only  non-linear 
conversions  necessary  are  simple  geometrical  ones  which  translate  from  focal 
plane  to  object  field  and  display  space.  The  advantages  of  working  in  the 
spatial  frequency  terms  will  become  clear  as  the  analysis  is  developed.  In 
the  following  development,  the  NVL  model  approach  will  be  used  precisely  but 
will  be  applied  in  the  foveal  lens  focal  plane  as  a function  of  linear 
spatial  frequency  (1/mm).  Parameters  will  be  covered  in  the  same  order  as 
they  are  in  the  NVL  Report  ( 4 ),  which  describes  the  model  in  detail. 

1.  MTF'S 

a.  Optical  MTF 

The  optical  MTF’s  consist  of  a diffraction  MTF  and  a Gaussian  MTF. 

(1)  Diffraction 

In  angular  terms,  the  diffraction  MTF  is  referenced  as  Equations  (9) 
and  (10)  of  the  NVL  report: 

Hopt(£x,0)  “ n [cos_1  A - A(1  - A2)1/2]  (D-l) 

where  A - XF^fx/L(6)  (D-2) 
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where  L(0)  is  the  equivalent  focal  length  which  changes  over  a 50/1  range  as 
object  field  angle  0 changes.  The  angle  6 is  the  absolute  angle  between  the 
point  of  interest  and  the  lens  optical  axis.  At  the  foveal  lens  image  plane 

f 

s - - y 

x L(0) 


(D-3) 


where  Sx  is  the  image  plane  spatial  frequency  and  f is  its  object  field 
angular  equivalent  measured  along  the  scan  line  projection  in  the  object 
field  (p  direction  on  Figure  D-l).  Solving  for  f in  Equation  (D-3)  and 
substituting  this  for  f^  in  Equation  (D-2).  p 


A - X F#  Sx  (D-4) 

Since  the  F number  of  our  lens  is  constant,  the  diffraction  MTF  is  no  longer  a 
function  of  object  field  angle.  Thus  we  may  write  H ^(S  ) which  indicates 
that  the  MTF  is  a function  of  the  independent  variab?iLS  xonly.  Note,  however, 
that  conversion  to  object  field  angular  spatial  frequency  is  very  simple 
because  focal  length  is  constant  over  small  angular  increments  and  may  be 
determined  from 


f = S L(0)  (D-5) 

p x 

where  p is  along  the  scan  line  projection  in  the  object  field 
likewise 

fw  = SyL(6)  (D-6) 

where  w is  normal  to  the  scan  direction  in  the  object  field 


(2)  Blur 


A similar  simplicity  exists  here.  The  MTF  equation  with  the  angular 
term  b of  Equation  (11)  of  Reference  ( 4 ) replaced  with  its  equivalent  is: 


- ~p[- S £*2]  (D-7> 

The  foveal  lens  inherently  has  a constant  spatial  blur  over  its  entire  focal  plane, 
so  that  the  sigma  (o)  of  Equation  (D-7)  Is  a constant.  Substituting  Equation  (D-5) 
into  (D-7)  we  see  the  blur  MTF  simplifies  to 

Hblur  ^Sx)  = CXp  2lf2  °2sx2]  (D-8) 

Thus  this  MTF  like  the  diffraction  MTF,  is  no  longer  a function  of  object  field 
angle  because  the  focal  length  variable  has  been  removed. 
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b.  Detection  MTF 

The  spatial  filter  MTF  of  the  detector  is  defined  as: 

Sin(nf  Ax) 

,JDet(fx‘0)  " irf^Ax  " Sinc(fxAx)  (D.9) 

It  is  also  complex  in  our  system  because  the  angular  projection  of  the  detector 
into  the  object  field  (Ax)  in  this  equation  varies  with  absolute  object  field 
angle  (0).  Since  the  detector  height  is  still  uniform  at  the  lens  focal  plane, 
shown  in  Figure  D-2  as  Ah,  Equation  (D-9)  can  be  restated  as: 

Sin(:rS  Ah  ) 

■wv  ■ —1  - <D-io> 

x x 

Again  the  MTF  becomes  independent  of  object  field  angle.  Note  from  Figure  D-2  that 
the  detector  height  (Ah^)  is  a function  of  detector  size(a),  detector  system  focal 

length  (1*D),  and  relay  focal  length  (Lc),  viz: 


Ah  - a Jfc  (D-ll) 

ld 

There  are  two  ways  to  arrive  at  the  required  L and  L based  upon  whether  the  scanner 

v#  l) 

or  detector  characteristics  are  known.  If  the  scanner  is  available,  horizontal  FOV 

(0  ) and  instantaneous  FOV  (a)  are  known.  Then  the  detector  focal  length  is 
6 


L 


D 


(D-12) 


On  the  forward  end  of  the  scanner  we  must  cover  the  foveal  lens  image  plane 
height  of  2b (6^)  which  subtends  the  angle  8 . Thus, 


Solving  for  Lc 


h(G  ) 

max 


LC  ’ 


h(e  ) 

max 

tan  0 /2 
8 


(D-13) 


(D-14) 


If  the  detector  characteristics  are  known,  the  focal  lengths  are  a function 
of  detector  size  (Ah)  projected  unto  the  image  plane  as  shown  in  Figure  D-2. 
Detector  size  Ah  can  be  computed  directly  from  either  the  on-axis  resolution 
required,  the  number  of  scan  lines  required  across  the  vertical  FOV,  or  bandwidth/ 
response  restrictions  and  frame  rate  requirements.  The  focal  lengths,  and  Lp, 

are  then  selected  to  make  the  detector  dimension  appear  as  the  required  Ah  at 

the  foveal  lens  focal  plane.  The  scanner  horizontal  FOV  given  by  scanner  mechanics 

must  cover  the  total  image  plane  height  (2h(0  )).  In  either  case  the.  detector 

max 
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Location) 

G»?t  07?)  ? 

Figure  D-2  Scanner  resolution  element  geometry 


MTF  becomes: 

,IDet(Sx)  = Sinc  f~lT 
x D 

Again  this  MTF  is  independent  of  object  field  angle. 


e.  Detector  Electronics  MTF 

It  is  in  this  MTF,  the  detector  electrical  response,  that  we  get  into 
real  trouble  trying  to  work  in  object  field  angular  space.  For  a conven- 
tional linear  optical  system,  a linear  detector  scan  velocity  converts 
into  a scaled  but  linear  angular  scan  in  the  object  field.  This  is  not 
true  in  our  system  as  was  shown  in  Figure  D-l.  A linear  scan  in  the  x 
direction  on  the  image  plane  results  in  angular  velocities  in  both  9X  and 
0y  directions  in  the  angular  object  field.  Both  of  these  angular  components 
are  non-linear  functions  of  both  x and  y position  on  the  image  plane.  Thus, 
converting  from  spatial  frequency  to  temporal  frequency  becomes  very  complex. 
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All  of  this  can  be  avoided  by  working  in  linear  spatial  plane  terms.  If 
the  scanner  has  an  angular  scan  velocity  then  the  linear  motion  of  the 
instantaneous  FOV  on  the  foveal  lens  image  is 


? 


V = 
x 


PI* 


The  conversion  to  temporal  frequency  (f)  is  therefore 


This  is  a constant  conversion  and  not  a function  of  time.  Therefore,  all 
electronic  MIT's  of  the  NVL  model  are  valid.  These  are 


“Elect  (f> 

H,(f> 

d.  Display 

The  RVS  display  is  the  inverse  of  the  foveal  lens,  which  results  in 
a conventional  linear  raster  generated  on  the  CRT.  The  CRT  has  a constant 
spot  size  and  the  expansion  optics  has  a constant  blur  at  the  object  focal 
plane.  Again  this  MTF,  if  derived  in  the  linear  spatial  plane,  will  not 
be  a function  of  object  angle.  If  the  optical  blur  and  CRT  spot  size  are 
combined  and  assumed  to  have  a Gaussian  MTF,  a composite  sigma  (04)  results 
and  the  MTF  is: 


WV  ' “xp  [-2"2(r<,d>  V ] 


where  r is  the  physical  ratio  of  format  sizes  ; viz 

„ “lens  image 
Display  crt 


By  contrast,  if  this  were  accomplished  in  the  object  angular  plane,  the  MTF 
would  be  much  more  complex,  viz 


2»2<r°d>Vl 

L(0)2M2  J 

where  M is  any  system  angular  magnification  from  object  field  to  the  viewer. 
Again  the  simplicity  is  obvious. 


HDisp<fx,0,M)  " exP 
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(D— 16) 

(D— 17) 


(D-18) 

(D-19) 


t 


e.  & f.  Stabilization  and  Eyeball 


The  remaining  two  MTF's  are  the  only  two  that  are  not  simplified  by 
working  in  linear  spatial  rather  than  angular  terms.  First,  stabilization 
tends  to  be  angular  input  to  the  system.  Using  the  MTF  from  the  NVL 
report : 


‘ '*>>  <-p£x2>  ®-20> 

Converting  to  the  foveal  lens  image  plane  results  in 

1L  (S  6)  - exp  [-PS  \(G)?]  (D— 21) 

JL«0  S X j X 


Similarly,  the  eye  views 


the  display  in  angular  terms. 


The  NVL  MTF  is 


Equation 


(D-22) 


must  be  converted  to 


the  foveal  lens  image  plane 


(D-22) 


(D-23) 


In  conclusion,  seven  MTF's  have  been  simplified  at  the  expense  of . two  that 
have  been  made  slightly  more  complex  by  the  conversion  to  linear  spatial 
frequency. 


2.  NEAT 

NEAT  can  be  used  as  defined  in  the  NVL  document  since  it  is  independent  of 

spatial  frequency,  therefore  no  changes  are  required.  However,  it  is  a function  of  F 

number  and  a discussion  is  in  order  concerning  which  is  the  correct  F number 

to  use  in  the  calculation.  It  is  not  simply  the  foveal  lens  F number  or 

objective  lens  F number  as  stated  in  the  NVL  document.  For  NEAT  calculations, 

the  F number  must  represent  the  actual  ray  cone  supplying  the  detector. 

In  general  our  lens  will  be  limited  to  a fixed  image  height  ^(®max) 

that  is  based  on  fabrication  considerations.  The  scanner  must  cover  this 
total  image  in  the  most  efficient  manner.  The  optical  geometry  of  the  lens 
•and  scanner  can  be  reduced  to  the  basic  arrangement  shown  in  Figure  D-3. 

On  this  figure  the  ray  cone  that  supplies  the  detector  is  <(>^.  Note  that 

a maximum  ray  cone  ($^)  that  the  detector  can  utilize  exists.  Obviously, 

the  best  design  would  have  4>y=4>^.  This  is  not  always  possible,  however, 

because  the  lens  may  be  required  to  operate  with  available  scanners.  Therefore, 
it  will  be  necessary  to  define  both  of  these  ray  cones  and  utilize  the  smallest 
for  NEAT  and  diffraction  MTF  calculations.  Before  processing,  note  that  ray 
cones  may  also  be  defined  by  F numbers  as  shown  in  Figure  D-4. 
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G*7ft  0773  4 


Figure  D-4  Ray  cone  parameters 


Since  anywhere  along  the  ray  bundle 


F#  = L/D 

$ - 2 tan_1(D/2L)  (D-24) 

and  F number  may  be  defined  by 


F 


it 


1 

2 tan  <f>/2 


(D-25) 


Therefore,  v:e  can  use  F number  instead  of  <J>  to  define  ray  cones.  By  the  nature 
of  collimntion,  the  ray  bundles  between  the  collimating  lens  and  the  detector 
will  be  of  constant  cross  section,  i.e.. 


1 

CM 

Pn 
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The  larger  of  these  two  F numbers  must  be  used  for  both  diffraction  and  NEAT 
calculations.  Now  note  that  F^  is  the  F number  of  the  foveal  lens  while  F^ 

is  that  generated  by  the  scanner  internal  aperture  stop  (Dg)  and  we  may  write: 
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(D— 28) 


For  any  particular  scanner,  Dg  will  be  known.  From  Equation  (D-27)  the  smallest 

F number  which  translates  into  the  largest  ray  cone  we  can  supply  the  detector 
with  is: 


(D— 29) 


which  is  also  obvious  from  inspection. 
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For  our  ONR  work,  a scanner  has  been  specified  while  for  the  USAF  effort  only 
a detector  geometry  is  available.  The  following  derivations  will  try  to  accoininodat 
both  arrangements.  If  the  scanner  is  specified,  Dg  and  will  be  known.  If  Lp 

is  not  given  it  can  be  calculated  from  the  scanner  instantaneous  FOV  (a)  and 
detector  size  (a)  in  the  same  plane.  Then 
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(D-30) 


Also,  if  the  scanner  is  specified,  its  maximum  scan  angle  (0g)  shown  in  Figure 

D-2,  is  known.  Since  the  scanner  must  scan  the  entire  image,  0g  must  cover 

2h(0  ).  Then 

max 
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(D-31) 


It  is  now  possible  to  compute  F^  from  Equation  (D-26). 
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where  F.,_  is  the  foveal  lens  F number.  Remember,  for  the  diffraction  MTF  and 
NEAT  calculations,  we  use  the  larger  of  F^  or  F^- 

If  a scanner  is  not  defined,  one.  must  be  theorized.  This  can  easily  be 
accomplished  by  procedures  outlined  in  Lloyd  (8  ).  Generally,  scanning  hardware 
(mirror,  drum,  etc.)  are  selected  to  accomplish  a given  scan  angle  0g  (number  of 

facets).  Collimation  lens  focal  length  (L^,)  can  then  be  computed  by  Equation  (D-31) 

From  the  required  scan  lines  across  the  maximum  image  height  h(0max),  the  apparent 

detector  size  (Ah)  at  the  image  can  be  calculated.  The  angular  instantaneous  FOV 
of  the  scanner  can  be  calculated  as  follows,  using  the  on-axis  geometry  of  Figure 
D-2: 
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Dg  can  then  be  defined  by  matching  the  lens  F number 
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3.  MRT  CALCULATIONS 

The  following  MRT  equation  modifications  are  required  so  that  the  computation 

may  be  performed  in  linear  spatial  frequency  terms.  First,  in  the  NVL  MRT  equation,  Ay 

must  be  replaced  by  the  apparent  detector  size  at  the  foveal  lens  image  plane, 

i.e.,  it  must  be  the  Ah  defined  on  Figure  D-2.  As  previously  demonstrated  in 

Equation  (D-ll) . 
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Also,  in  the  MET  equation,  it  is  best  to  compute  the  Q integral  in  terms  of 
temporal  frequency.  This  eliminates  the  velocity  term  in  the  MET  equation 
and  makes  the.  Q integral  easier  to  compute.  The  Q integral  is  therefore 


Q(f.,)  -/"i 


(D— 37) 


Of  these  terms,  only  H^,  the  transfer  function  for  a rectangular  bar  of  width  w, 

has  not  been  defined.  This  transfer  function  is  in  linear  rather  than  angular 
dimensions,  i.e.. 
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The  MRT  equation  written  to  show  the  dependency  of  two  variables  is 
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This  equation  results  in  an  MRT  very  weakly  dependent  on  0.  To  obtain  the  MPT 
for  any  field  angle  0 we  convert  the  spatial  frequency  term  into  an  angular 

frequency  term  by  using  Equation  (D-5)  containing  the  focal  length  function: 

fy  " SXL<6>  (D-5) 

Note  this  will  be  the  angular  spatial  frequency  in  the  scan  direction  (target 

bars  normal  to  the  scan  direction).  It  could  be  related  to  f and  f but 

X y 

this  docs  not  appear  to  be  required  at  this  point. 
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